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FOREWORD 

This User's Manual has been prepared to provide the engineer with the 
information required to run the coupled mode version of the Normal Modes 
Rotor Aeroelastic Analysis Computer Program. The manual provides a full 
set of instructions for nanning the program, including calculation of 
blade modes, calculacions of variable induced velocity distribution and the 
calculation of the time history of the response for either a single blade or a 
complete rotor with am airframe (the latter with constant inflow) , The 
analysis used is discussed in the technical report 1 Reference a. 


f 


- I 

k 




29 REV 0 


11 

PAOE 


I 




k 


i 


t 





Sikorsky RIrcraft 


U 


ofviwQN uMmio A»»CMArr oo wp o watiqn 



REPORT NO. SER-50910 


TABLE OF CONTENTS 


FOREWORD 

DISCUSSION 


What the Program Does . 
How This is Accomplished 
Natural Frequencies 
Rotor Steady State 
Helicopter Trim . . 
Transient Mode . . , 
Simplified Flow Chart 
What the User Must Do 


SECTION 1 
SECTION 2 
SECTION 3 
SECTION U 
SECTION 5 
SECTION 6 
SECTION 7 
SECTION 8 
SECTION 9 
SECTION 10 
SECTION 11 
REFERENCES 


BASIC DECK SET UP 

DETAILED DESCRIPTION OF DECK SET UP 

NUMERICAL LIST OF LOADER INPUT ITEMS 

LOGICAL LIST OF LOADER INPUT ITEMS 

SUMMARY OF PROGRAM CONTROL OPTIONS 

DETAILED DESCRIPTION OF LOADER INFUT ITEMS . . . . 

COMMON ERRORS AND PITFALLS 

ERROR MESSAGES AND DIACMOSTICS 

GENERAL SUGGESTIONS FOR RUNNING PROGRAM 

ACCURACY 

OUTPUT PRINTOUT 


Page 

ii 

1 

1 

2 

2 

2 

2 

3 

h 

5 

6 
7 

15 

21 

27 

28 
52 

52 

5J» 

5U 

55 

66 


jLU 


PAO£ 


$A KCV 0 
/ 



Sikorsky PIrcraft 


u 

OfVWiaN UNTTVO AIMCAA^ BCW»nW»T10M 


REPORT NO. SER-50910 


DISCUSSION 


This guide will answer these three questions. What is the program capable of 
doing? How, in general terms, is this accomplished. What must the user do to 
obtain the desired results. The emphasis is on how. The why is described in 
References a, b, c, d, and e. 

WHAT THE PROGRAM DOES 


The function of the program may be logically considered in five parts. 

The first part calcula.es the coupled natural frequencies and mode shapes from 
the blade physical properties and rotor speed. It is possible to stop at 
this point and run successive cases varying rotor speed or other blade 
physical properties to investigate blade frequency characteristics. 

The second part calculates rotor steady state, air loaning, blade bending 
and rotor generated forces and moments for a specified flight condition and 
control setting. Ihe program may be stopped at this point and multiple 
cases run to investigate parametric variations in control positions, blade 
p'nysical properties, or fli^t conditions. It also writes data files used in 
part three. 

The third part is the calculation of the veuriable inflow distribution across 
the rotor aisc. This part requires the blade response for a. steady state 
revolution to provide the necessary input data. The program can run either 
with constant inflow assumed or with the calculation of variable inflow. This 
part also generates a data file of induced velocities that is read, by tha 
coupled mode dynamic analysis section if induced velocities are to be included 
in the blade dynamic response. 

A full description of this section may be found in Appendix C, including the 
stand alone option. 

The fourth part is the trimming mode. In trim, equilibrium is established 
for all three orthogonal forces and all three orthogonal moments. It is 
possible to stop at this point and run successive cases. 

The fifth part is the calculation of the aeroelastic transient during a 
gust encounter. Again successive cases may be run. It is not possible to run 
variable induced velocity during this mode. 
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HOW IS THIS ACCOMPLISHED 


DISCUSSION (cont’d) 


Natural Freouencies and Mode Shapes 

The equations for the calculation of coupled flatwise and chordwise natural 
frequencies and mode shapes are described in Reference -b..' The equations 
for torsional frequencies and mode shapes are displayed in Reference 4. 

There are two hinge conditions permitted in the flatwise and the chordwise 
direction; articulated or nonarticulated. In torsion provision for either a 
root spring or rigid root attachment has been provided. 

Rotor Steady State 

The rotor blade modal amplitudes and first derivatives are given at the start 
of the first revolution of the rotor in the program. These starting values 
may be assumed zero, be left over from a previous case, or be loaded in via 
cards. The equations of motion described in Reference c are then employed to 
calculate the azimuthal time history of the blade motions. Three types of blade 
airfoil section data may be selected by the user; spar data - univariant function 
of angle of attack, bivariant function of angle of attack and Mach number, 
quadre.varient Sanction of the angle of attack and its first and second time 
derivatives, and Mach number. 

A reference blade is assumed at zero azimuth with the prescrj.bed starting 
values. The equations of motion of this blade are integrated around the 
azimuth from zero to 360 degrees. The proced’ore used is described in Reference 
dv At this point the values of modal eimplitudes and first derivatives at zero 
degrees azimuth and 360 degrees azimuth are compared. If the modal amplitudes 
and first derivatives are within a user specified tolerance the calculation 
proceeds to the integration of the rotor forces and moments and the calculation 
of blade bending moments and other output data. If convergence is not obtained, 
successive revolutions are performed until convergence is obtained or the maxi- 
mum number of revolutions, as specified by the user, is exceeded. If the 
number of revolutions is exceeded without convergence, the calculation proceeds 
after printing a warning that the blade mode shapes did not converge. 

Helicopter Trim 

By the use of partial derivatives and a Newton-Raphson iteration new values 
of lateral and longitudinal cyclic pitch, main and tail rotor collective pitch, 
and body roll attitude cure obtained and the Y force, Z force, and the moments 
about the X, Y, and Z socis are zeroed out. No attempt is made to zero the X 
force because of its -nonlinear chaoracteristics. However, some control of this 
force is provided by the provision for loading in the an^e between the rotor 
shaft and the relative wind. This angle remains essentially unchanged during 
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DISCUSSION (cont'd) 

successive iterations. The sequence of events is as follows: The rotor blade 

motions are brought into steady state for initial values of cyclic and 
collective pitch. The rotor forces and moments and blade flapping are cal- 
culated, also the body fox’ces and moments are calculated. The rotor Z force 
and helicopter gross weight are summed and compared with zero. If this in 
addition to blade flapping is also zero within the tolerances specified by 
the user, the tail rotor collective pitch and body roll attitude are 
calculated to zero the yawing moment and the side force. The calculation 
then proceeds to a new case or the transient section. If the rotor Z force 
does not balance the gross weight or the flapping is not zero and the maximum 
number of trim iterations has not been exceeded, the cyclic and collective 
pitches are adjusted and the rotor is brought to the steady state associated 
with the new control positions. This procedure continues until trim 
convergence is obtained or the maximum number of iterations is exceeded. 

If the case does not converge, after printing an appropriate warning, the 
calculation proceeds as in the converged situation. 

Newton's method is used to predict new control positions at the end of each 
iteration. This method uses the nine partial derivatives asscci4te4 with 
rotor Z force and flapping with respect to cyclic and collective pitch. The 
procedure is described in detail in the technical manual. The partial deri- 
vatives can be calculated internally by the program, either by paraitetric 
permutations or by linear theory. If known, provision has been made for 
loading them in via cards. This, if possible, will significantly reduce 
the compute time. 

Transient Mode 


Tne last section is the transient calculation and gust response. After 
helicopter trim has been attained the helicopter may be subjective to a 
vertical g'ust with a gust front of any arbitrary shape. Penetration of 
both the body and the rotor disc may be simulated. During the transient 
mode the response of each blade is calculated individually euid the 
reference blade may be assumed at any azimuth position at the time the gust 
front and the rotor disc first come into contact. The gust can penetrate 
the disc at any angle and at any velocity. The gust front velocity is in- 
dependent of the helicopter forward velocity. The actual velocity of gust 
penetration is calculated internally. At each azimv.th position up to a 
maximum, specified by the user, blade stresses, deflections, and body motions 
are printed out. Vfhen the maximum azimuth angle (number of revolutions) is 
attained the program either terminates or proceeds to a new case as requested 
by the user. 
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Siuaiaary of What the Prograa Does and Hov It Does It 

A simple flow chart may be used to demonstrate the five basic sections. 



There are two basic modes of program operation. The first is the sane as 

Deck. IlUl which calculates the aeroelastic response of a helicopter rotor 

and support structure in a wind txmnel. This method of operation is described • 

in Reference d. The second is the inclusion of variable induced velocity, | 

The equat-tons for the calculation of variable induced velocity emd a detailed ^ 

description of the computer program may be found in Appendix C. The variable 

induced Valbcitv is completely coupled euad no additional input cards are j 

required. The procedure is as follows. The procedure is divided into three 

steps. Step 1 first caLLculates the coupled natural frequencies and mode 

shapes if not loaded via card input then calculates rotor and/or helicopter 

steady state utilizing the coupled blade dynamic equations as given in Refer* 

ence c . Data, files of bdeic input data and output from this case are written 

ftor use in steps 2 and 3. Step 2 calcxilates the actual inflow across the 

rotor disc utilizing the data files written during step 1. The equations 

used are given in the references of Appendix C. When the induced velocities 

have been calculated, they are written on a data file to be read by the 

following step. Step 3 reloads the basic blade/helicopter data that was 

written during step 1 and reads the variable Induced velocities that were 

written at the end of step 2. The coupled mode rotor blade response is then 

calculated and the rotor /helicopter retrinued again using the coupled mode 
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DISCUSSION 


dynamic equations as in Step 1. No f’urther data need be loaded except 
**JCL^* and control cards. Provision for modifying the data between steps has 
been provided by two individual cpJ.1s to loader for each step. These options 
are normally used for debugging arid the program could be modified to eliminate 
these double calls if desired. There is no need for the two calls to loader 
in step 1, but the second call was not bypassed to be consistent. 

What the User Must Do 


The user has four basic functions to perform. 

1. Set up proper system control cards. 

2. Load in all data in the proper format. 

3 . Be sure that what is asked for is really what is wanted. 
h. Be able to recognize errors if they occur auid make the 

appropriate corrections. 
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SECTION 1 


mASIC DECK SET UP 


1 . Rijin card 

2. JCL for step 1 

3. Ir-put data for step 1. 

U. -1 99 -1- 

5. Additional input data 

6. -1 99 -1- , .V 

7 Curve Data (If card option used; 

s! Mode Shapes (If loaded via cards) 

9. JCL for step 2. 

10. Control Card 

12*. ChaSea^to Basic Data Generated by step 1 if desired 

13. -1 99 -1. 

lU. JCL for step 3 
15 -1 99 2. 

16. Changes to input data read in step 1 

17. -1 99 -1- ^ ^ 

18. Curve data (if loaded via cards) 

19. JCL for termination 
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SECTION 2 DETAIL DESCRIPTION OF DECK SETUP 

1. RUN CARD 

Defined by system used. 

2. JCL FOR STEP 1. 


The JCL for all three steps and termination may be found In Appendix A. 
Slight differences between systems may exist. 

3. INPUT DATA FOR STEP 1 


Card input data for the case being run in loader format. Muca of the input 
data may be loaded in shotgun form. (Unlike the curve data which is loaded 
in block fora no firm order to input need be adhered to.) The input para- 
meters that are likely to change between cases are of this nature. Each 
item is assigned a location number and may be loaded singularly if desired. 
The location-parameter allocation is supplied in Sections 3 and U. A 
detailed descriptior. of the items is supplied in Section 5* The loader 
format is as follows: 


mu-Ji - II ’•8-J 


JEi- 


JL 


I 2 3 


67 


l£ 19 


30 31 


42 43 




JL 


54 ss 


..£ 9 . 



T - to end case 

N number of parameters to be loaded in on a card. Maximum of 

five. 

L Location of PI. P 2 is loaded into L+1. P3 inti L+2, etc. 

Right adjusted. 

PI Numerical value of PI anywhere between columns 7 and I 8 . 

P2 Numerical value of P2 anywhere between columns 19 and 30. 

P 3 » P^» P5« etc. is right adjusted and has no decimal point. 

PI - P^ must have decimal points. 

N May be any number from 1 to 3. 


If L is not specified, PI is assumed to follow the last location specified 
by the previous card. 

A saapls of the shotgun input data appear*' in the Appendix. 
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SECTION 2 


DETAIL DESCRIPTION OF DECK SETUP (cont»d) 


U. FIRST END OF DATA CARD 

•♦1 99 + 1. In loader format indicates end of data - begin calculation. If 
location 522 contains a 0., Induced velocities are not to be calculated. 

'*1. indict^tes subsequent cases to follow. . indicates tt mlnation after 
current case has been completed. If location 522 contains a +1. data files 
are to be written for step 2. 

5 . ADDITIONAL DATA 


This is only Included for consistency with the rest of the s'- up. There /ould 
normally not be any cards loaded. 

6. SECOND END OF CARD 
Same comment as Item U. 

7. CURVE DATA 


Curve data loaded in CLOAD format or from data file 22. If location 135 
contains A ^ 1. data •'“lie 22 is used. This must be preloaded with the 
section data. If location 135 contains A >1- data is loaded via cards. 

A significant waount of inp\it is not normally changed between cases during a 
specific computer run. The blade aerodynamic section data and the body 
aerodynamic data are good examples. The following is a complete list of 
the curve data that may, depending on requested options, be needed. These 
data are loaded by means of subrout.' ne "CLOAD". 

The following symbols are used. 

Local blade section aerodynamic angle of attack - degrees 

N Local blade section Mach No. 

A Function of Op (see Reference a) 

B Function ofHp (see Reference a) 

Body pitch angle of attack positive nose up - degrees 

a Body yaw angle of attack positive clockwise - degrees 


Wing angle of attack positive pitch up - degrees 
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SEC non 2 

! DETAIL DESCRIPTION 

OF DECK SETUP (cont’d) 

CUP.\^ ID 

DESIGNATION OUTPUT VARIABLE 

UI!IVARIANT BIVARIANT TRIVARIANT 


Load If Spar Data 

Is Reauestad 


CLSP 

CL 



CDSP 

CD 

“r 


C^6P 

CM 

<»R 


Load 

If Steady Two Dimensional or Unsteady Section Data 

la Requested 

CLDAT 

CL 

®R 

M 

CLDDAT 

CD 

®R 

M 

Ct®AT 

CM 

®R 

M 


Additional Data If Unsteady Secti'cn Data Option Is 

Requested 

CNUNST 

CN 

“R 

A ] 

cMasT 

CM 

®R 

A ] 

ASTCN 

CL at stall 

M 


/u3TCM 

CM at stall 

M 


RCNiM 

C!lg^ at Mach/CNg^ at Kach » 0 

M 


CN'REGx 

Lower Boundary a cf CN 
Transition Region 

M 


CHFBG2 

Upper Bo’ondwyoj^ of CN 
Transition Region 

M 


CFr^EGl 

Lower Boundary® ^ of CM 
Transition Region 

M 


CMREG2 

Upper Boundary® r of CM 
Transition Region; 

H 



(A (Xill explanation of the mechanics cf the unsteady aerodynamics 
option Is presented In Reference f). 
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SECTION 2 

DETAIL DESCRIPTION OF DECK SETUP (eont'd) 


Load if Body Force Option is Requested 

CURVE ID 

VARIABLE 

UNIVARIANT 

ENORM 

Body normal force coefficient 

“'f 

BLONG 

Body longitudinal force 
coefficient 

° F 

BSIDE 

Body side force coefficient 


BPIT 

Body pitching moment coefficient 

«F 

BYAW 

Body yawing moment coefficient 

®F 

qTAIL 

Tail dynamic pressure ratio 

« F 

DWTAIL 

Tail downvash 

Generalized nondimens ional main 
rotor slipstream penetration 
distance 

NF«T 

Normal force coefficient 
horizont6Ll tail 

Horizontal tail angle of 
attack degrees positive pitch 
up 

NFVT 

NormaJ. force coefficient 
vertical tail 

Vertical tail angle of attack 
degrees positive yaw clockwise 


Load if WinK Option is 

Selected 

WNOIF 

Coefficient of wing normal force 

dy 

WNGCP 

Coefficient of wing chordwiae 
force 

®W 

WNGPM 

Coefficient of wing pitching 
moment 

®W 

WNGROL 

Coefficient of wing rolling 
moment 



Load if Gust Option is 

Selected 

GUST 

Value of vertical gust velocity 
in ft /sec 

distfluice of gust penetration 
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SECTION 2 Detail Description of Deck Set Up (cont'd) 

The following three curves have been added to facilitate 
the use of the coupled variable induced velocity option. They need not 
be loaded if induced velocities are not to bo calcxilated. 
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SECTION 2 DETAIL DESCRIPTION 07 DECK SETUP (con<-»d) 

Load it Varia'ble Indn'’^ed Option is Used 
ACLZ ALPHA for Zero Cj^ vs Mach No. 

ACLM ALPHA at Maxlmiin Cj, vs Mach No. 

ASLP Lift Curve Slope CL/RAD vs Mach No. 

END 

While a specific order of the curves is not necessary for operation of the 
program, a logical order, such as that given above, is desirable for good 
organization, curves pertaining to options not requested need not be loaded. 
Curves loaded but not referenced do not interfere with program operations. 

CLOAD Format 

Curve data is input using subroutine CLOAD. This load curves in a 
standard, easily handled format, which is defined as follows; 
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SECTION 

2 

DETAIL 

DESCRIPTION 

OF DECK SETUP (cont’d) 

1. 

Alpheunaneric entries are indicated by A. 

2. 

Decimal entries are 

indicated by F must have decimal point. 

Entry 

Type 

Card Number Column Number 

Comment 

ID 

A 

1 

1-6 

Curve identification 

Title of ciurve as defined in data 

TITLE 

A 

2 

1-72 

statement of calling progrsun 

X TITIJC 

A 

3 

I- 2 U 

X axis title 

Y TITLE 

A 

3 

25-1*8 

Y axis title 

IJPTS 

A 

u 

1-1* 

Coded literally, NPTS actually 
pii.iched in columns 1-U 

N 

F 

k 

8-10 

Number of X-Y pairs corresponding 
to one bi variant value 

ALPHA 

A 

h 

1 I +-19 

Bivariant title 

ALPHA 

F 

k 

21-30 

Bi variant value 

BETA 

A 

u 

3 U -39 

Trivariant title 

BETA 

F 

k 

1 * 1-50 

Trivariant value 

X,Y etc 

. F 

5 

1-10,11-20 

X-Y pairs (N of them) 

etc . 

• 



etc. thru 80 

« 

• 

END 


last card 

1-3 

Coded literally, actually punch 


END in columns 1-^3 • 

Only if this is last cuz^s. 


After the X-Y pairs have been «ntered« car-’.a k and 5 Are continued for all 
parsuneter changes, hivariant varying fastest. If there is another curve, 
vrlte curve identification and continue. *END* is coded only vhen all 
curves are entered. 

A sample list out of lift coefficient in univariant, hivariant and tri- 
variant from is displayed in Appendix. A.' . 
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SECTION 2 DETAIL DESCRIPTION OF DECK SETUP (cont'd) 

8. MODE SHAPES 

If location 101 contains a ~i. node shapes are loaded via cards previously 
punched from the mode shape module. If location 101 contains a +1. mode 
shapes are calculated. 

9 . JCL FOR STEP 2. 

See Appendix A for details. 

10. CONTROL CARD FOR INDUCED VELOCITY 

Normally no cards loaded. See Appendix C for detailed information. 

11. END OF DATA 

-1 99 -1. is loaded. 

12. CHANGES TO DATA LOADED FROM STEP 1. 

Normally no cards loaded. Can be used for debugging options. See Appendix C. 

13 . END OF DATA 
Load -1 99 ~1. 

lU. JCL FOR STEP 3. 

See Appendix A for sample JCL, 

15 . STEP 3. CONTROL CARD 

To execute step 3 the following card must be loaded -1 99 2. (Loader Format), 

16 . CHANGES TO DIPUT LO A DED IN STEP 1. 

Normally no cards loa&ed. Can be used for debugging options. 

17. END OF DAT A 

Load -1 99 •• (Loader Format). 

18 . CURVE DATA (If loaded via cards) 

Same comment aul Itei ’f. 
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SECTION 2 DETAIL DESCRIPTION OF DECK SETUP (cont'd) 

19 . JCL FOR TERMINA' ^ ION 
See Appendix A* 


Mode Shapes for Step 3 . 

Unlike the curve data which must be reloaded, the mode shapes have been 
written on a data file during step 1. and need not be loaded. 

SECTION 3 


A list of shotg\in input items in ascending numerical order follows. This 
list is most useful in setting up original (first time) runs or for deter- 
mining what a particular item is from its location number. 
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LOADER 

LOCATION INPUT ITEM UNITS 


1 

2 

3 

k 

5(15) 

20(15) 

35(15) 

50(15) 

65(15) 

80(15) 

95 

96 

97 

98 

99 
100 
101 
102 
103 
lOl* 

105 

106-124 

125 

126 

127-128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 
1»*5 

146 

147 
l4a 
149 
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OMECiA-R (TIP SPEED) FT/SEC 

RADIUS FT 

OFFSET RATIO ND/ RADIUS 

YOUNG'S MODULUS PSI 

SEGMENT LENGTH ND/RADIUS 

SEGMENT MASS (INCLUDING COUNTERWEIGHT) SLUGS 

FLATWISE INERTIA IN»*4 

CHORDWISE INERTIA IN*»4 

TORSIONAL INERTIA LBS-3EC*»2 

TORSIONAL STIFFNESS (GJ; LBS-1NS»*2 


NUMBER OF BLADE MODES USED 
NOT USED 
NOT USED 

DEBUGGING OPTION DYNAMICS 
=-l. TO TERMINATE JOB 

TORSIONAL ROOT SPRING IN-LBS 

=^.CALC MODES, READ MODES 

NON ZERO TO PUNCH I-K)DES 

=1. TO CALCULATE MODE SHAPES ONLY 

FLAWISS =0. ,ARTIC =1.,N0NARTIC 

NOT USED 

NOT USED 


EDGEWISE =0 . ,ARTIC =1.,N0NARTIC 

PRINT HELICOPTra; FORCES IP' +1'.' ' ' 

NOT PART OF LOADER VECTOR 

1«1AXI:TJM TRANSIENT AZIMUTHAL ANGLE 

NUMBER OF BLADES 

TIP LOSS 

INCRE!*ENTAL DRAG COEFFICIENT 
=1.F0R NO AERODYNAMICS ON SEGMENT ONE 
NUMBER OP SEGMENTS (MUST BE 15.) 

-1. DATA ON CARDS +1. DATA ON^FILB ' 

LINEAR WIST IN AERODYNAMIC MODULE 

GRAVjTY 

PRE-CONING 

PRE-LAGGING 

TAIi(DELTA3) 

TAN (ALPHA!) 

NOT USED 
LAG DAMPER 

CHORDWISE STRUCTURAL DAMPING 

NOT PART OF LOADER VECTOR 

=1.T0 PRINT ROTOR STEADY STATE TIME HISTORY 

TOTAL NONLINEAR TWIST IN AERODYNAMIC MODULE 

NOT PART OF LOADER VECTOR 

»1.T0 RUN HELICOPTER TRANSIENT 


DEGS 


DEGS 

FT/SEC*»2 

DEGS 

DEGS 


LBS-SEC-FT/SEC 

ifl). 


4 




DEGS 


-1 
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150(15) 

CHORD 

FEET 

165 ( 15 ) 

HOIILINEAR TWIST IN AERODYNAMIC MODULE 

DECS 

180 ( 15 ) 

I/C FLATWISE 

1N»*3 

195(15) 

I/C CHORDWISE 

IN»»3 

210(15) 

NOT USED 


225(15) 

BUILT IN rWIST USED IN NATURAL FREQUENCY CALC. 

DEGREES 

2i+0(l5) 

DIST QUARTERCHORD IS FORWARD OF ELASTIC AXIS 

nd/radius 

255(30) 

NOT USED 


285 ( 15 ) 

COUNTERWEIGHT MASS 

SLUGS 

300(15) 

NOT USED 


315(15) 

DIST CNTRWT IS FORWARD OF ELASTIC AXIS 

ND/RADIUS 

330 ( 15 ) 

BUILT IN TWIST FIRST DERIV 

RAD/N.D.SPAH 

31 * 5 ( 15 ) 

BUILT IN TWIST SECOND DERIV 

RAD/N.D.SPAN**2 

360 ( 15 ) 

DIST SPAR CENTROID FORWARD OF ELASTIC AXIS 

HD/RADIUS 

375(60) 

NOT USED 


1*35(15) 

DIST BLADE CG IS FORWARD OF ELASTIC AXIS 

ND/RADIUS 

450 ( 15 ) 

HOT USED 


1*65 

AIR DENSITY 

SLUGS/PT^‘»3 

1*66 

SPEED OF SOUND 

FT/SEC 

1*67 

SHAPE TOLERANCE 

ND 

1*68 

NUMBER OF SHAPE TRIALS 


1 * 69 - 1*98 

NOT PART OF LOADER VECTOR 


1*99 

WARNING READ DESCRIPTION- NO. OF HELICOPTER TRIM 

ITERATIONS 

500-506 

NUT PART OF LOADER VECTOR 


507 

ALPHA .ANGLE MADE BY SHAFT WITH RELATIVE WIND 

DECS 

508 

NOT PART OF LOADER VECTOR 


509 

T0LEH.M1CE ON Z FORCE 

LBS 

510 

TOLERANCE ON A13 FLAPPING 

DEG 

511 

TOLERANCE ON BIS FLAPPING 

DEG 

512 

NOT PART OF LOADER VECTOR 


513 

AIS 

DECS 

51 I* 

BIS 

DECS 

515 

THETA 75 ; 

DECS 

516 

INFLOW RATIO, IJVMBDA 

ND/OMEGAH 

jl7 

HELICOPTER FORWARD VELOCITY 

KNOX’S 

518 

AZItvIUT: INTERVAL (DYNAMICS) 

DECS 

519 

AZIMUTH INTERVAL (PRINT AND INTEGRATE) 

DECS 

52C 

»1.F0R SONSY!^«!ETRICAL AIRFOIL 


521 

NOT PART OF LOADER VECOTR 


522 

»1.T0 USE COUPLED VARIABLE INDUCED VELOCITIES 


523 

»1.T0 PRINT COMPLETE TRANSIENT OF ROTOR 


5"1» 

NOT USED 


525 

»1.T0 PRINT NAT FREQUENCIES AND MODES DEBUG 


526 

-l.TO PRINT AERODYNAMICS DEBUG 


527-532 

NOT PART OP LOADER VECTOR 


533 

TRANSIENT PRINT CONTROL. SEE TEXT 


53 U 

CURVE NUMBER USED BY SEGME2JT 1 


535 

CURVE hL*MBER USED BY SEGMENT 2 


536 

CURVE NUMBER USED BY SEGMENT 3 


537 

CURVE NUMBER USED BY SEGMENT U 
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538 

CURVE NUMBER USED BY SEGMENT 5 


539 

CURVE NUt'SER USED BY SEGf-ENT 6 


5J4O 

CURVE NUl-IBER USED BY SECrffiNT 7 


5 ia 

CURVE NUMBER USED BY SEGIENT 8 


5U2 

CURVE NU^3ER USED BY SEG!<ENT 9 


5U3 

CURVE NUMBER USED BY SEGMENT 10 


541* 

CURVE NUMBER USED BY SEGt-SNT 11 


5^*5 

CURVE NU>®ER USED BY SEGIENT 12 


546 

CURVE NUMBER USED BY SEGf-ENT 13 


547 

CURVE NUMBER USED BY SEGt-ENT l4 


548 

CURVE NUMBER USED BY SEGfENT 15 


549 

=1.T0 CALCULATE BODY FORCES 


550 

=1. COUPLE FUSELAGE MOTIONS TO ROTOR BLADE 


551-552 

NOT PART OF LOADER VECTOR 


553 

BUTTLINE f4AIN ROTOR 

INS 

554 

DOWNWASH FACTOR , FUSELAGE 

ND 

555 

FUSELAGE STATION , CG 

INS 

556 

WATERLINE , CG 

INS 

557 

FUSELAGE STATION , BODY DATA REF. 

i.*r. 

558 

WATERLINE , BODY DATA REF. 

INS 

559 

P'USELAGE STATION , HORIZ TAIL AERO CENTER 

INS 

560 

WATERLINE , HORIZ TAIL AERO CENTER 

INS 

561 

FUSELAGE STATION , MAIN ROTOR 

INS 

562 

WATERLINE , MAIN ROTOR 

INS 

563 

MAIN ROTOR LONGITUDINAL SIIAPT TILT 

DEG5 

564 

LATERAL DISPLACEMENT, TAIL C OF PRESSURE 

FEET 

565 

HORIZONTAL TAIL AREA 

FEET»»2 

566 

HORIZONTAiL TAIL INCIDENCE ANGLE 

DECS 

5f'l 

VERTICAL TAIL AREA 

FEET*»2 

568 

FUSELAGE STATION , VERT TAIL AERO CENTER 

INS 

569 

WATERLINE , VERT TAIL AERO CENTER 

INS 

570 

WINGSPAN 

FEET 

571 

FUSELAGE STATION , WING AERODYNAMIC CENTER 

INS 

572 

WATERLINE , WING AERODYNAMIC CENTER 

INS 

573 

WING AREA 

FEET»»2 

574 

WING INCIDENCE ANGLE 

DECS 

575 

ROLL DAMPING TERM , WING 

ND 

576 

FUSELAGE STATION , TAIL ROTOR HUB 

INS 

577 

WATERLINE , TAIL ROTOR HUB 

INS 

578 

BUTTLINE , TAIL ROTOR HUB 

INS 

579 

TAIL ROTOR CANT ANGLE 

DECS 

580 

NUMBER OF BLADES , TAIL ROTOR 


581 

TAIL ROTOR CHORD 

FEET 

582 

TAIL ROTOR RADIUS 

FEET 

583 

TAIL ROTOR ANGULAR VELOCITY 

RAD/8EC 

584 

NOT PART OF LOADER VECTOR 


585 

TAN DEL3 TAIL ROTOR (NOT USElf AT PKESEOT) 

ND 

536 

TIP LOSS , TAIL ROTOR 


587 

TAIL ROTOR BLADE INERTIA 

SLUG-FT»»2 
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588 

589 

590 

591 

592 

593 
59»* 

595 

596 

597 

598-604 

605(10) 

615(10) 

625-799 

'iOO 

801-812 

313 

814 

815 

816-318 

819 

820 
82] 

822-824 

825 

626 

827 

828-636 

837 

838 
339 

840 

841 

842 

843 

844 

345 

846 

847-350 

8=^1 

852 

353 

854 

855 

856 

857 

858 

859 

860 


DOWKWASH FACTOR , TAIL ROTOR 

AIRCRAFT GROSS WEIQiT 

POLAR ^10MEriT OF IHERTIA ABOUT X AXIS 

POLAR t-Dt^IIT OF INERTIA ABOUT Y AXIS 

POLAR MOMENT OF INERTIA ABOUT Z AXIS 

X-Z PRODUCT OF INERl'IA 

FRACTION OF GROSS WEIGHT CARRIED BY WING 

WING Mi-'AN AERODYNAMIC CHORD 


LATERAL SHAFT TILT 

MASS MCl-ENT , TAIL ROTOR BLADE 

NOT USED 

INITIAL VALUE QW 
INITIAL VALUE QW<^ 

NOT USED 

=1. CALCULATE HARMONICS OF BLADE BENDING ANGLE 
NOT PART OF LOADER VECTOR 
PARTIAL OF Z P'ORCE WITH RESPECT TO AlS 

PARTI.AL OF Z FORCE WITH RESPECT TO BIS 

PARTIAL OF Z P'ORCE WITH RESPECT TO TH75 


NOT PART OF LOADER VECTOR 
PARTIAL OF AlS FLAPPING WITH 
PARTIAL OF AlS FLAPPING WITH 
PARTIAL OF AlS FI>APPING WITH 
NOT PART OF LOADER VECTOR 
PARTIAI. OF BIS FLAPPING WITH 
PARTIAL OF BIS FLAPPING WITH 
PARTIAL OF BIS PXAPPING WITH 
NOT PART OF LOADER VECTOR 
1. FOR GUST 

VELOCITY OF GUST FRONT 
ANGLE MADE BY PERP WITH X DIRECTION 
AZIMUTH POSITION OF REF BLADE AT INCIDENCE 
AZ INCREMENT FOR PRINT OF AZIP!UTHAL STRESS 


RSPCT-AIS FEATHERING 
RSPCT BIS PEATHERING 
RSPCT TH75 

RSPCT AlS FEATHERING 
RSPCT BIS FEATHERING 
RSPCT ra75 


LBS 

SLUG-FT»»2 
SLUG-FT«*2 
SLUG FT**2 
SLUG-FT««2 

FEET 

DEGS 

SLUG-FT 


FT/SE' 

DEGS 

DEGS 

DEGS 


=l.TO PRINT AZ iVERO INFO IN TRANSIENT 

=l.TO PRINT BLADE MOTIONS IN TRANSIENT 

VALUE OF LAG DAMPER OF REFERENCE BLADE DURING TRANSIENT 

TH75/NP DEGS 

TWl/NF DEGS 

NOT USED 

TOTAI. NL uJER OF BODY FREEDOMS 
NOT USED 

NUMBER OF RETAINED MODES 
ROTOR HEAD TRANSLATIONAL MODE NO -X 
ROTOR HEAD TRATJSLATIONAL f40DE NO -Y 
ROTOR HEAD TRANSLATIONAL MODE NO -Z 
ROTOR HEAD ROTATIONAL MODE NO - ABOUT X 
ROTOR HEAD ROTATIONAL FWDE HO - ABOUT Y 
ROTOR HEAD ROTATIONAL MODE NO - ABOUT Z 
STRUCTURAL DAMPING 
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861-869 not used 

870 l.TO INPUT M, K 

2 . TO INPUT M, K, G 

3 . TO INPUT MD, KD, G 

871 O.TO INPUT STIFFNESS, K 

l.TO INPUT FLEnaiLITY 

872 l.TO PRINT EIGENVECTORS AND VALUES 

873 l.TO PRI.NT ALL MATRICES 

87 ^+ l.TO PRINT BODY VELS, ACCS 

875 l.TO PRINT MODES 

876 l.TO PRINT BODY DISPLACEt-ENTS 

900 l.TO USE FLEXIBLE BODY MODES 
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Section U 

A list of shotgun input items grouped by function follows. This list is 
useful for making changes to existing complete deck setups, or if the program 
is used to calculate natural frequencies and mode shapes only. 
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IWPUT PJ\RAMErrERS III LOGICAL ORDER 
NATURAL FREQI'EIICIES AND i\K)DE SIIAPES 


LOCATION IIJPUT ITEM 


UNITS 


1 TIP SPEED 

2 RADIUS 

3 OFP’SET RATIO 

U YOUNG'S MODULUS 

5 ( 15 ) SEGMEirr LENGTH 

20 ( 15 ) SEGMEIiT M;\SS (INCLUDING COUNTERWEIGHT) 

35 ( 15 ) FLATWISE INERTIA 

50U5) CHORDWISS INERTIA 

65 ( 15 ) TORSIONAL INERTIA 

80 ( 15 ) TORSIONAL STIFFNESS 

95 BLADE f-KDES USED 

100 T0R3ICNAI. ROOT SPRING 

102 - 1 . TO PUi.'CH ^DDES 

103 =1. TO CiVLCULATJ .’DDE SHAPES ONLY 

10** FLArVISE =0. ,ARTIC =1.,N0NARTIC 

125 EDGF^ISE =0., ARTIC =l.,NONARTIC 

525 =1. TO PRINT NAT FREQUENCIES AIID MODES DEBUG 

101 *-l. TO READ MODES =+l. TO CALCULATE KfiODES 


FT/SEC 

FT 

ND/RADIUS 

PSI 

ND/RADIUS 

SLUGS 

IN*»1* 

LBS-SEC**2 

LBS-INS»*2 

IN-LBS 


PROGRAM CONTROL OPTIONS 

U 99 WARNING READ DESCRIPTION- NO. OF HELICOPTER TRIM ITERATIONS 

800 *1. CALCULATE HARMONICS OF BLADE ROOT FLAPPING 

135 ’-1. CURVE DATA ON CARDS *+l. ON DATA FILE 

5^*9 =1. TO C/vLCULATE BODY FORCES 

550 »1. COUPLE F-USELAGE J-OTIONS TO ROTOR BLADE 

59*‘ FRACTION OF GROSS WEIGHT CARRIED BY WING 

837 =1. FOR GUST 

ll»9 »1. TO RUN HELICOPTER TRANSIENT 

99 ■-!. TO TERMINATE JOB AND CONTROL STEP WHEN USING VARIABLE INDUCED 

VELOCITIES 

522 -1. TO USE VARIABLE INDUCED VELOCITIES 

101 »-l. TO READ MODES =♦!. TO CALCULATE MODES 

102 ■!. TO PUNCH MODES 

lOU FLATWISE »0.. ARTIC -1., HONARTIC 

125 EDGEWISE »0.. ARTIC »1.. NOKARTIC 

103 »1. TO CALCULATE MODE SHAPES ONLY t ♦1. IN 525 ALSO 
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PRINT OPTIONS 


lue 

=1. TO PRINT ROTOR STEADY STATE TIME HISTORY 


519 

AZIf-f’rH INTERVill. (PRINT ,1ND INTE'IRATE) 

DECS 


= 1. 'J'' I HIN ■’ Ci;.‘-'i"ijET E 1 h/ui J Ir.N L UE R'jlOR 


535 

TPA'l’IEir: PRINT CON PROD, SEE USER GUIDE 


0]i 

AN r.iCPlMENT RiR PRINT n- A '-’l Mr" HA'. S'i’''ESS 

DF.rs 


= 1. :o PRINT AZ AERO INFO IN TRANSIFDJT 



BIJVL’E PARAME'fERS 


u’» 

NUMSKR OF SEGMENTS (MUST BE 15) 


2 

RAolIUS 

FT 

3 

C-rS RATIO 

:ID/RAD1US 

120 

N'lMJER Or BIjVDES 


5(15) 

s- ::t-:::t li-::: oru 

KD/RA.OIUS 

^0(15) 

SEGIENT M\SS (INCLUDING COUNTERVEIGHT) 

SLUGS 

150(15) 

CHORD 

FECT 

1 Ji -- 

Ail AiiiAU.)DY.» AMI s MODULE 

'T'£GS 

13. 5) 

N>.N.INTj\R T-;IST in A.^'RODYNAMLC MODULE 

DEGS 

llT 

T' '.'al no;iline>\r ivist in aerodynamic j^dule 

DEGS 

132 

PRE-CjNINQ 

DECS 

139 

PIuE-LAGGING 

DEGS 

l'*3 

LAI DAMPER 

LBS-SEC-FT/RAD 

61»U 

VA!,i-K or LAG DAMPER Of REFERENCE BLADE DURING TRANSIENT 

154 

CHTl-.r.'A'ISE STRIJ.OrURAL DAMI’ING 

ND 

l-'to 

TAN (DELIA 3) 


lul 

TLNiALIHA 1) 


1.00(15) 

I/C EIjVIVISE 

IIl»*3 

195(15) 

I/C CHORDWISE 

IN**3 

•:10(15) 

DIST DIMKTERCKOHO IS FO'WARD OF ELASTIC AXIS 

ND/RADIUS 

315(15) 

Dir.T CMRaV IS lORVARD OF ELASTIC AXIS 

ND/RADIUS 

330(15) 

BUILT IN TWIST, FIRST DERIV. 

RAD/ND SPAM 

3‘.5(15) 

BUILT IN TWIST, SEC. DERIV. 

RAD/ND SPAN»*2 

360 ( 15 ) 

DIST SPAR CFINTROID FORWARD OF ELASTIC AXIS 

ND/RADIUS 

^35(15) 

DICT BLADE CG IS FORWARD OF ELASTIC AXIS 

SD/RADIUf* 


NUMBER Or SHAPE TRIALS 


467 

SHAPE TOLERANCE 

ND 

518 

AZItrjTH INTERVAL (DYNAMICS) 

DECS 

519 

AZirrj’TH INl'ERVAL (PRINT AND INTEGRATE) 

DEGS 

605 ( 10 ) 

INiriAL VALUE QW 


615 ( 10 ) 

INITIAL VALUE QW» 
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AERODYNAI-HC OPTIONS 


131 TIP LOSS 

133 =1. FOR NO AERODYNAMICS ON SEGMENT ONE 

132 IKCREMFNTAL DRAG COEFFICIENT 

520 =1. FOR NONSYMMETRICAL AIRFOIL 


534 

CURVE 

NTJI-IBER 

USED 

BY 

SEGt^ENT 

1 

535 

CURVE 

NUMBER 

USED 

BY 

SEGMEOT 

2 

536 

CURVE 

NUMBER 

USED 

BY 

SEGMENT 

3 

537 

CURVE 

NUMBER 

USED 

BY 

SEGMENT 

4 

538 

CURVE 

NUMBER 

USED 

BY 

SEGMENT 

5 

539 

CURVE 

NUJ.I3ER 

USED 

BY 

SEGMEirr 

6 

5k0 

CURVE 

NUISER 

USED 

BY 

SEGIffiNT 

7 

5hl 

CIWE 

NUI-BER 

USED 

BY 

SEGMENT 

8 

5k2 

CURVE 

NUI-ffiER 

USED 

BY 

SEGNIENT 

9 

543 

CURVE 

NUIffiER 

USED 

BY 

SEGMENT 

10 

544 

CURVE 

N'JMBER 

USED 

BY 

SEGMENT 

11 

545 

CURVE 

NUI©ER 

USED 

BY 

SEGMENT 

12 

546 

CURVE 

NUNSSR 

USED 

BY 

SEGMENT 

13 

547 

CURVE 

NUMBER 

USED 

BY 

SEGt^ENT 

14 

548 

CURVE 

NUMBER 

USED 

BY 

SEGMENT 

15 


ROTOR AND FLIGHT CONDITIONS 


1 

TIP SPEED 

PT/SEC 

513 

AIS 

DEGS 

5l4 

BIS 

DECS 

515 

THETA 75 

DEGS 

507 

ANGLE MADE BY SHAFT WITH RELATIVE WIND 

DEGS 

5l6 

INFLOW RATIO 

ND/OMEGAR 

517 

HELICOPTER P’ORWARD VELOCITY 

KNOTS 

466 

SPEED OF SOUND 

PT/SEC 

465 

AIR DENSITY 

SLUGS/PT»»3 

137 

GRAVITY 

PT/SEC**2 
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FUSELAGE PARAI-IOTEBS 


5U9 =^1. TO C;\LCUIJ^TE BODY FORCES 

550 = 1 . COUPLE FUSEIJVGE MOTIONS TO ROTOR BLADE 

563 LONGITUDINAL SHAFT TILT, MAIN ROTOR 

596 LATERAL SHAFT TILT, MAIN ROTOR 

589 AIRCRA.'-T GROSS WEIGHT 

590 POLAR KO'-IENT OF INERTIA ABOUT X AXIS 

591 POLAR MCjMENT of INERTIA ABOUT Y AXIS 

592 POLAR !-' 0 MENT OF INERTIA ABOUT Z AXIS 

593 X-Z PRODUCT OF INERTIA 


DECS 

DECS 

LBS 

SLUG-PT »*2 

SLUG-FT *»2 

SLUG-FT »*2 

SLUG~FT **2 


( 




553 BUTTLINS ^ 4 ,AIN ROTOR 

555 FUSEIA^GE STATION, CG 

556 WATERLINE , CG 

557 FUSEIAGS STATION, BODY DATA REF. 

558 WATERLIiJE , BODY DATA REF. 

561 FUSEL,\GE STATION, MAIN ROTOR 

562 WATERLINE , MAIN ROTOR 

55U DOWirWASH FACTOR, FUSELAGE 


INS 

INS 

INS 

INS 

INS 

INS 

INS 

ND 


59U FRACTION OF GROSS WEIGHT CARRIED BY WING 

573 WING AREA 

570 WINGSPAll 

595 WING tiEAN AERODYNAMIC CHORD 

WING INCIDE.NCE A:^GLE 

575 roll da:^ping term, wing 

571 FUSELAGE STATION , WING AF:R 0 DYKAMIC CENTER 

572 WATERLINE , WING AERODYNAMIC CENTER 


FEBT *»2 

FEET 

FEET 

DEGS 

INS 

INS 


565 HORIZONTAL TAIL AREA FEET »»2 

566 HORIZONTAL TAIL INCIDENCE ANGLE DEGS 

561 « LATERAL DISPLACEMENT, TAIL OF PRESSURE FEET 

559 FUSELAGE STATION, HORIZ TAIL AERO CENTER INS 

560 WATERLINE , HORIZ TAIL AERO CENTER INS 


567 VERTICAL TAIL AREA FEET »»2 

568 FUSELAGE STATION, VERT TAIL AERO CENTER INS 

569 WATERLINE , VERT TAIL AERO CENTER INS 
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576 

FUSEUGE STATION, TAIL ROTOR HUB 

INS 

577 

VVATERLINE , TAIL ROTOR HUB 

INS 

578 

BUTTLIKE , TAIL ROTOR HUB 

INS 

580 

NUMBER OF BLADES, TAIL ROTOR 


581 

TAIL ROTOR CHORD 

FEET 

582 

TAIL ROTOR RADIUS 

FEET 

583 

TAIL ROTOR ANGULAR VELOCITY 

RAD/SEC 

586 

TIP LOSS , TAIL ROTOR 


587 

TAIL ROTOR BLADE INERTIA 

SLUG-FT«*2 

597 

MASS IviOMENT, TAIL ROTOR BLADE 

SLUG- FT 

579 

TAIL ROTOR CANT ANGLE 

DEGS 

588 

DOWNWASH FACTOR, TAIL ROTOR 




liELICOFTER STEADY STATE 


U99 

WARNING READ DESCRIPTION-NO. 

OF HELICOPTER TRIM ITERATIONS 

589 

AIRCRAFT GROSS WEIGHT 



LBS 

511 

TOLERAI1CE ON BIS FTIAPPING 



DEG 

510 

TOLERANCE ON AIS FLAPPING 



DEG 

509 

TOLERANCE ON Z FORCE 



LBS 

819 

PARTIAL OF alG FLAPPING WITH 

RSPCT 

.AlS FSLATHSRING 


820 

PARTIAL OF alS FLAPPING WITH 

RSPCT 

BIS FEATHERING 


821 

PARTIAL OF alS FLAPPING WITH 

RSPCT 

TH75 


825 

PARTIAL OF bis FLAPPING WITH 

RSPCT 

AlS FEATHERING 


826 

PARTIAL OF bis FLAPPING WITH 

RSPCT 

BIS FEATHERING 


827 

PARTIAL OF b LS FLAPPING WITH 

RSPCT 

TH75 


813 

DZ/DAIS 




8l4 

DZ/DBIS 




815 

DZ/DTH75 






GUST PARAMETERS 


837 

=1. FOR GUST 


838 

VELOCITY OF GUST FRONT 

FT/SEC 

839 

ANGLE MADE BY PERP WITH X DIRECTION 

DEGS 

8U0 

AZIMUTH POSITION OF REF BLADE AT INCIDENCE 

DEGS 

129 

MAXIMUM TRANSIENT AZIMUPHAL ANGLE 

DEGS 


DEBUG OPTIONS 


525 

»1. TO PRINT NAT FREQUENCIES AND MODES DEBUG 


526 

«•!. TO PRINT AERODYNAl'JiCS DEBUG 


98 

» TO PRINT ROTOR BLADE DYNAMICS DEBUGGING 



2. TO PRIHT A. S, and T MATRIXES SEE REFERMCE AT EACH 
DELTA FSI PLUS ABOVE ITEMS 

3. PRINT COMPONENT PARTS MAKING UP ELEMENTS OF S, AND T. 
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SECTION 5 SltW^ARY CF PROGRAM CONTROL OPTIONS 


The following description defines in general terras what must he done to run 
the various options of the progran. 

Location 101 controls the manner in which the node shapes are made available 
to the program. A +1.0 indicates that mode shapes are to be calculated by the 
program. A -1.0 indicates that mode shapes are to he loaded in via punched 
cards. A 0.0 indicates that mode shapes previously loaded or calculated are 
to be used. If this location is zero for the first case a default option of 
+1.0 is used. This location is reset to zero after each case. Therefore, 
wnenever new mode shapes are to be used this location must be ±1.0. Mode 
shapes and frequencies as punched by the program are in the proper form and 
may be read in directly. 

Rotor Blade Steady State 

The limit on the number of rotor revolutions for the establishment of rotor 
steady state is loaded into location * 168 . Locations 605 - 62 U have been pro- 
vided for the starting values of the modal ampliti’des and first derivatives. 

If these initial guesses are close to the actual .steady state values the 
remount of compute time will be reduced. If the values are urduiown, zero should 
be loaded. 

Rotor Trim Iteration 

The maximum number of trails that the program is allowed to converg on rotor 
Z force and flapping is loaded into location ^^99. The procedure used is 
described in the technical report. Reference a, and the data needed is de- 
scribed in Section U of this manual. If a zero is loaded into this location 
no iteration is performed. 

Body Forces and Helicopter Trim Iteration 

This also includes the uing and the tail assembly aerodynamic forces. This 
option is exercised by loading a +1.0 into location 5-‘9. If this option 5s 
requested body physical constants, see Section U, and body aerodynamic data, 
see Item 10. of Section 2, must also be loaded. The addition of the body 
forces has no effect on the rotor iteration. But tail rotor collective and 
body roll attitude are adjusted to balance yawing moment and side force. 

If a +1.0 is loaded into location 550 the fuselage motions are couple to the 
blade dynamics. This is normally done if body forces are calculated. The 
loading of +1.0 into locations 5l*9 and 550 switch the program from a rotor 
simulation to a helicopter simulation. 
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Transient Running 

If a transient response is to be calculated this may be done by loading a 
+1.0 into location 1U9. It is also necessary to load a maximtmi value of 
azimuthal angle inxo location 129. This angle is in degrees. One 
revolution is S60 degrees, two revolutions, 720, etc. This angle does not 
include any trim iterations. 

Gust Response 

Penetration of a vertical gust may be simulated by loading a +1.0 into 
location 837. Other necessary gust parameters are given in Section It. 

Section 6 - Detailed Description of Loader Input Items by Location Number 

1. Tip Speed 

Main rotcr tip speed expressed in feet/sec. 

2. Radius 


Main rotor radius measured from center of rotation expressed in feet. 

3 . Offset Ratio 


Distance flapping and lagging hinges are displaced from the center of 
rotation nondimensionalized by the radius. 

U. Young's Modulus 

2 

Young's modulus expressed in Ibs/ln . 

5 . Segment Lengths (15 ) 

For eterodynamic and elastic considerations the rotor blade is considered 
to be composed of fifteen segments defined from root to tip, each having 
constant physical and aerodynamic properties. Although these segments 
can be any length desired, it is good practice to make the segments shorter 
outboard since the aerodynamic forces become nonlinear and are of greater 
relative magnitude. It is also common practice to make the last (Segment 
i.^} equal to 1. > tip loss; see location 131. The segment lengths are 
nondimensionalized by the radius. Thus, it is essential that the sum of 
the segment lengths plus the offset ratio equal 1.0 i .0001. 
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20. SepT,ent Maas (l5) 

The mass of each segaent expressed in slugs. This includes the mass of the 
counterweights. 

35" Flatwise Area Moment (15) 

Flatwise area monent of inertia about the elastic axis of each segment 
expressed in inches^. 

50. Chordvise Area Moment (15 ) 

Chordwise area moment of inertia about the elastic axis of each segment 
expressed in inches^. 

65. Torsional Inertia ill) 

Torsional mass moment of inertia about the elastic axis of each segment 
expressed in (Slug-rt^)/ft. (Mote this is also expressible as Ib-Sec^.) 

80. Torsional Stiffness (l3 ) 

Torsional stiffness of each seffient expressed in lbs-in2. 

95 • Number of Blade Modes (Frequencies ) 

The desired number of blade modes (including rigid body, flap, and lag 
modes, when applicable). Modes are included in order of increasing 
frequency. Up to 7 modes may be included. 

98. Debugging 

Nonzero in this location activates a debugging printout of the equations of 
motion. 

99* End of Case/Job 

A minus word count indicates the end of a case. A minus entry indicates 
the end of the run after the completion of the last case. If location 522 
contains a -*’1. (variable induced velocity option) , this location should 
contain a +1. for Step 1, or a +2. for Step 3. 

100. Torsional Root Spring 

Value of torsional spring In Inch-lbs between blade and rotor head. If 
rigidly attached 1.0 x 10^^ Is recommended. This value may be loaded In 
as 1.0E20 (E type format), provided it is right adjusted in the parameter 
field of a normal loader card. 
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101. Mode Shape Calculation Control 

0.0 Do not calculate or read mode shapes and frequencies. 

+1.0 Calculate frequencies and mode shapes. 

-1.0 Load mode shapes and frequencies from cards. 

This location must contain +1., if location 103 contains +1. for each case. 
This location is automatically set to zero after mode shapes have been 
calculated or loaded except when locationK)3 has a +1. to stop calculation 
after obtaining modes. Thus the second and subsequent cases will use the 
mode shapes and frequencies used by the first case. If at any time new 
mode shapes are desired, this may be accomplished by loading at ^1.0 into 
this location. Again upon completion of the case, this control will be 
reset to zero. 


102 . Punch Modes 

A +1.0 in this location activates a routine that punches the mode shapes 
and derivatives in a format compatible with this program. 


103. Mode Shape Only Option 

A +1.0 in this location indicates that the program is to stop after the 
calculation of frequencies and mode shapes. In this case locations 101 and 525 
must also contain a +1. for each case. 


lOU. Flatwise Rotor Hinge Condition 

Load 0.0 for articulated, +1.0 for nonarticxilated. 


125. Chordvlse Hinge Condition 

0.0 articulated chordwise +1.0 nonarticulated chordwise. 
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126 . PebuFiging of Equations of Helicopter Motion 

A +1. in this location prints all forces acting on helicopter. Normally 
activated. 

127-128 


Not part of loader vector. 

129 . Maxlr.un Transient Azimuth Angle 

The caximum value azimuth angle, PSI, may attain while in trwslent mode. 
For example, for two revolutions read in 720. 

130. Number of Blades Main Rotor 

Limited to six if transient option is exercised. Any number may be used 
if transient not required. 


131. Tip Loss, B 


Tip loss defines the length of the blade effective in producing lift. 
Normally tip loss is 0.97* However, in order to take into account all 
possibilities, provision is made for any tip loss. The CL and the CM of 
the last (Tip) segment are adjusted by the following formulas: 


C 




(DX^^ -1. + B)/DX^^ 





+ b)/dx^^ 


It may be seen 


from these formulas that if DX 


15 


1- B; 



and 


C„ are zero. This condition will result in the best aerodynamic 
15 

accuracy. Unfortunately, however, if B is closer to 1. than .98, this will 
result in a of .02 or less. This may result in some inaccuracy in 

the calculation of the natural frequencies and mode shapes. For this 
reason delta-x*s less than .02 are not recommended. In this case, make 
DX^j» .02 and let the program adjust CL and CM. Note DX^^j. must be greater 

than l.*-B. When using variable inflow tip loss is generally set to 1.0. 
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132. Delta Drag 

An incremental drag coefficient added to all aerodynamic drag coefficients 
obtained from the blade section properties. £;apirically determined for a 
specific airfoil for the 0012 .002 is normally used. 

133. Root Loss Option 

A + 1 . loaded into this location omits all aerodynamic forces on the first 
' Root ) segment . 

13U. Number of Segments 

Ttiis is not an option. The program requires that 15 . he loaded into this 
location. 

136. Linear TVlst 

The total amount of built-in linear aerodynamic twist from center of 
rotation to blade tip expressed in degrees. Normally negative. If twist 
is nonlinear omit. If blade is not twisted, 0.0 must be loaded. 

137 . Gravity 

2 

Gravity expressed in fi/sec . 

138 . Preconing 

Amount of precor.ing for nonarticulated rotors expressed in degrees . 
Positive up. 

139 . Prelagging 

Amount of prelagging for nonarticulated rotors expressed in degrees. 
Positive in direction opposite to rotation. 

140. Tangent of Delta-3 

Tan of Delta-3. A positive value reduces the blade pitch aa a linear 
function of blsuie root angle (flapping). 

lUl. Teuigent of Alpha 1 

Tan of Alpha 1. A positive value increases pitch as blade lags. 
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II4 3 . Lae Dampin g 

Amount of lead/lag root end damping expressed in ft-lbs-sec/rad. This 
does not affect chordwise frequencies in this analysis. 

llU. Chordvise Structural Dampin g 

Amount of chordwise structural damping of nondimensional form. Values 
in the order of +0.02 to +O.OU are normally used. This does not affect 
the chordwise frequencies. 


1I45. 

Not part of loader vector. 
li*6. Print Time History 

A +1.0 loaded into this location results in a complete printout of the 
steady state (or the last attempt to achieve steady state) blade modal 
amplitudes and first and second derivatives. 

148 . 

Not part of loader vector at present. Skip or load in zero. 

149 . Calculate Transient 

A +1.0 loaded into this location indicates that after helicopter trim has 
been accomplished or the maximum number of iterations has been exhausted 
a transient response considering each blade individually is desired. 


150 . Chord Main Rotor (l?) 


The chord of each segment root to tip, expressed in feet. 


165. Nonlinear Twist in Aerodynamic Module (19) 

If the twist is nonlinear the twist value of each segment is loaded in 
root to tip expressed in degrees. This twist should be referenced to the 
.T5R radius station. The twist is used in the aerodynamic module. Twist 
values for structural variations are loaded in locations 225 - 239« 
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180. I/C Flatwise (15) 

Hie flatwise blade area moment of inertia divided by the distance from 
the elastic axis to the extreme flatwise structural fiber expressed in 
inches^, loaded from root to tip. These values are used only to obtain 
stress from blade moments calculated by the program. 

1Q5. I/C Chordwlse (15) 

The chordwise blade area moment of inertia divided by the distance from the 
elastic axis to the extreme chordwise structural fiber expressed in 
inches^. Loaded for each segment root to tip. These values are used 
only to obtain stress from blade moments calculated by the program. 

225. Built-In Twist — Nonlinear 

The built-in twist values for each segment is loaded in root to tip and 
expressed in degrees. This twist should be referenced to the .75R radiua 
station. These values are used in the natural frequency calculations. 

These values need not be loaded if they are the same as the aerodynamic 
nonlinear twist. 

2J+O. ETA C/t* (15) 

The cho-dwise distance from the quarter chord to the elastic axis, 
nondimensionalized by radius. Loaded for each segment root to tip. 

Positive for the quarter chord ahead of the elastic axis. 

285. Counterweight Mass (15) 

Mass of the counten/etght of each segment expressed in slugs, loaded root 
to tip-.. 

300. Not used 

315 • Counterweight I>;catlon (15) 

Chordwise distance from counterweight to elastic axis , nondimensionalized 
by radius loaded for each segment. Root to tip. Positive for each 
counterweight ahead of elastic axis. 

330. Built-in Twist First Derivative (15). Nondimensional 

The first derivative with respect to radius distance (span) of the built-in 
twist. Local value of each station expressed as radians/nondlihensional 
span. Loaded root to tip. Must be loaded if either aerodynamic or built 
in twist is nonlinear. 
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3 )j 5. Butlt-in TwiBt Second Derivative (13 ) 

Second derivative of "built-in twist expressed as radlans/nondimenslonal 
span squared. Loaded only if twist is nonlinear. Must be loaded if either 
built in or aerodynanic twist is nonlinear. 

360 . Location of Spar Centroid (15 ) 

Chordwlse distance from spar centroid to elastic axis, for each segment, 
nondimensionalized by radius. Loaded root to tip. Positive for spar 
centroid oiiead of elastic axis. 


375* Not used 

390. Not used 

U05. Not used 

U20. Not used 

1*35- Location of Blade C.G. (13) 

Chordwise distance from blade center of gravity to elastic a.xis, 
iiondiaeasionalized by radius. Positive for C.G. forward of elastic 
axis. 

U50. Not used 

U65 . Air Density 

3 

Air density in slugs/ft . 

Ii66 . Speed of Sound 
Speed of sound in ft/sec. 

U67 . Mode Shape Tolerance 


Convergence tolerance on flapping, lee'ing, and modal amplitudes and 
first derivatives. Normally a vadue of 0.001 is used. A tolerance five 
times larger is automatically used when lead and lead prime are being 
tested for convergence. 

U68. Shape Trials 

The number of rotor revolutions for mode shape convergence. If thin 
number of rctor revolutions does not result in repeated values (within 
the specified tolerance) at azimuthal positions of 0.0 and 360, the 
>■ calculation proceeds after printing a warning that the blade node shapes 

I did not converge. 

I 
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U69-U98. 

Hot part of loader vector. 

’* 99 . Number of Trim Iterations 

The maximum number of iteration trials alloved to converge on Z forco 

and a, anl b, if articulated, or Z force and first harmonics of first 
Is Is 

flatwise bending amplitude, if nonarticulated. If the paraieters are not 
trimmed within the specified toler.ances after th.^s n’lmber of attempts a 
warning so stating is printed and the calculation continues. 

If this number is positive the program calculates the partial derivatives 
that ore used in predicting new control positions. Since this requires 
a substantial amount of computer time, these derivatives may be loaded 
directly into the program. A negative number of trim iterations indicates 
that the partials are to be loaded. If the number of trials is negative, 
but all partials are zero, the partials are calculated by linear theory. 
(See Reference a. ) 

If this option is used tail rotor parameters must be loaded. This does 
not apply to the INJIVAC IIO8/IIIO. 

500-506. 

Not part of loader vector. 


507. Angle Made by Shaft with Relative Wind 

If the transient option or variable induced velocity is requested, it is 
necessary to load in the inclination of the rotor shaft relative to the 
direction of flight, 'fhis angle is positive for a shaft tilted aft. 

This exercises some control on trim in the X direction. The units are 
expressed in degrees. This should not be confused with the physical 
shaft longitudinal tilt with respect to body (location 563). 
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508 , 


Aot part of loader vector. 

The coruiient not part of loader vector indicates that the particular 
L .:cation(s ) are not us»_d for input in the particular pro^/ram. Girce 
t:;ore are Govora] versions of the normal moles pro^^ram with inter- 
•j:i'i:uj,^ral)lo or r*oarly interchaiirtcabl o inputs, it is not , 70 od practice 
IS ■'.oiiv'y any projcram an. i this Ic'cation to a new inpuL parameter. 

I’uis location can be sicipped or loaded with zero with the r?mt of the 
vvar.e information. 
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509-511. Trim Tolerances 

The following locations are used to input the convergence tolerances on 
rotor forces and moments. 

509 Z force lbs 

510 a flapping deg. 

J. s 

511 flapping deg. 

513- 

Longii^udinal cyclic pitch degrees, negative harmonic series, AIS cosine 
coefi icier t. 

51U. BIS 

Lateral cyclic pitch degrees, negative harmonic series, BIS sine coefficient 


515. Theta-75 

Collective pitch at .75R degrees. 


517. Helicopter Forvard Velocity, V 
Forward velocity in knots. 


518. Delta-PSI Dynamics 

Azimuthal integration interval in degrees used in the calculation of 
blade dynamics. Constraint 1: The value must be evenly divisible into 

360. Constraint 2 : (Only applicable if treinsient option is exercised). 

Value must be evenly divisible into the angular spacing of the blades. 
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Although there are no other physical restrictions on this parameter, 
particular care must he exercised in its selection since the computing 
time is inversely proportional to its size. Too large a value will 
result in a mathematical instability, while too small a value will 
result in long compute times • A handle on the optimum value may he found 
by dividing the highest blade I'requency in cycles/rev (flatwise, 
chordwise, or torsional) into 36. The quotient obtained may then he 
rounded up to the nearest number that satisfies the defined constraints. 

In practice, this is often 5 deg. For 5-bladed rotors, this value would 
be 1+ deg. It should be pointed out that if many nodes are requested for 
stiff rotor blades, this parameter could be 1 deg. or less. This results 
in exceptionally long compute times. Normally a delta PSI less than 
2 degrees indicates that too many modes have been requested. 

519 . Delta-PSI Force Integration 

This Delta PSI is used in the calculation of the rotor forces in the 
printout. This Delta PSI must be selected such that (Delta PSI integrate)/ 
(Delta PSI dynamic) is an integer. Also the value must not be less than 
five degrees. 

520. Nonaymnetric Airfoil 

If this location is 0.0 the program will assume the airfoil is symmetrical 
and expect the blade section data to he defined only for positive angles 
of attack. If a + 1.0 is loaded into this location the blade is assumed 
to be nonsymmetric and the program requires blade section data for both 
positive and negative angles of attack of .-the spar and steady data, 
fhe unsteady data is assumed symmetrical. 

521. Not part of loader vector. 

522. Coupled Variable Induced Velocities to be Calculated 

A + 1.0 in this location results in the writing of files to be used in the 
calculation of vsuriable induced velocities. Note Alpha(s) (l«catioB 
307) must also be loaded. This is done in Step one. No further action 
is needed in Step three. 

523. Print Full-Blade Time History 

A + 1.0 in this location prints a complete time history of blade mod a l 
amplitudes and derivatives at every azimuthal station, diuring both the 
trimming and transient modes. Hiis results in a great deal of output 
and is primarily used for debugging. 

32U. Not used. 
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525 • Debugging of Natural Frequency ajid Mode ShaT?e3 

A +1.0 in this location results in a printout of the iteration used to 
determine natural frequencies. Also must he +1, if mode shape only option 
is requested. 

526 . Debug Aerodynamics 

A +1.0 in this location results in a printout of the components that make 
up the aerodynamic forces and moments. 

527 - 532. 

Not part of loader vector. 

533. Print Azimuth Information 

'A -1.0 in this location vill result in a printout of blade aeroelastic 
and aerodynataic information at each azimuth position as defined by location 
519» for both the trimming and the transient helicopter modes. A +1.0 
suppresses the printout in the trimming mode. There are two further 
controls available to reduce printout in the tremsient mode; see locations 
8^1 and 842. A zero suppresses all azimuth printout. 

534 - 548. Section Data Option (15 ) 

In the calculation of C. » and there are three selections of airfoil 

Li/ N 

data tables (if tables are loaded) available for each segment. This is 
controlled by the nmnber indicated as follovs: 

-*-1.0 Use univariant data table f (angle of attack) 

+2.0 Use bivariant data table f (<*, Mach number (M)) 

+3.0 Use unsteady data table ffa,a,(i) * 

0.0 or not loaded defaults to +2.0 (bivarlant data). Loaded for each 
segment, any mixture is permitted. For loading in of airfoil data, see 
Section 2. 

5U9. Calculate Body Forces 

A +1.0 in this location indicates that body forces are to be calculated. 

For a helicopter gust study this must be used. If 0. is loaded body forces 
are not cedculated and body input items need not be loaded. 

550. Couple Rotor and Body 

A +1.0 in this location indicates that the rotor and the body are to be 
coupled. For a helicopter gust study this should be used. Implementation 
of this option feeds hub motions into the blade aerodynamics. 
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551 - 552. 

Not part of loader vector. 

553. Buttline Main Rotor 

Lateral distance from body plane of symmetry to location of center of 
rotor hub. Expressed in inches, positive to ri^t. 

55^ . Fuselage Dovnvash Factor 

Ratio of average dovnvash at fuselage to momentum value at rotor. 
(Nominal value 0.6). 

555. Fuselage Station C.G . 

Fuselage station of body center of gravity expressed in Inches. 

556 . Waterline C.G . 

Waterline of body center of gravity expressed in Inches. 

557 . Fuselage Station. Body Data Ref . ^ 

Fuselage station vhere body data are defined expressed in inches. 
Fuselage stations dimensions increase frcsn nose to tail. 

558 . Waterline, Bo dy Data Ref . 

Waterline vhere body data are defined expressed in inches. Waterline 
dimensions increase from bottom to top of aircraft. 

559 . Fuselage Station of Horizontal Tail Center of Pressure 
Expressed in inches. 

560 . Waterline of Horizontal Tail Center of Pressure 
Expressed in feet. 

561 . Fuselage Station Center of Main Rotor Hub 
Expressed in inches. 

562 . Waterline Center of Main Rotor Hub 
Expressed in inches. 
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563. Main Rotor Shaft Tilt Relative to Body 
Expressed In degrees* positive forvard 

56U. Lateral Displacement Horizontal Tail Center of Pressure 
Expressed in feet* positive ri^t. 

565. Horizontal Tail Area 

Square feet. For no horizontal tail use zero. 

566. Horizontal Tail Incidence Angle 
Expressed in degrees. Positive leading edge 19. 

567. Vertical Tail Area 

Expressed in square feet. For no vertical tall use zero. 

568. Fuselage Station of Vertical Tail Center of Pressure 
Expressed in Inches. 

569. Waterline of Vertical Tail Center of Pressure 
Expressed in inches. 

570. Wing Span 

Tip to tip e 3 q>ressed in feet. 

571. Fuselage Station of Wing Center of Pressure 
Expressed in inches. 

572. Waterline of Wing Center of Pressure 
Expressed in inches. 

573. Wing Area 

Square feet (total for hoth vlngs). 

37U. Wing Incidence Angle 

Expressed in degrees. Positive leading edge up. 
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575* Wing Roll Dacping Factor 

NondimensioneLl roll damping contribution of the wing. It may be evaluated 
in terms of wing geometry using Figure 11 of HACA Report IO98, 01*, 
alternatively from DATCOM. Nominal value is -0.35* 

576. Fuselage Station of Tail Rotor Hub 
Expressed in inches. 

577. Waterline of Tail Rotor Hub 
Expressed in inches. 

578. Buttline of Tail Rotor Hub 
Expressed in inches. Positive to right. 

579. Tail Rotor Cant Angle 

Expressed in degrees. Positive tilt to ri^t. 

580 . Number of Blades of Tail Rotor 

Since linear theory is iised for the tail rotor calculation, any msnber 
of blades may be used. 

581. Chord of Tail Rotor 

Chord is assumed constant, expressed in feet. 

582. Radius of Tail Rotor 
Expressed in feet. 

583. Rotational Velocity Tail Rotor 
Expressed in radians/sec. 

58U. 


Not part of loader vector. 

585. 


( 





Hot part of loader vector. 
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586 . Tj-p LosSiTall Rotor 
Nominal value 0.97« 

58T. Second Maas Moment of Inertia of Tall Rotor Blade 

2 

Expressed in slijg-ft . 

588. Downvash Factor »Tall Rotor 

Ratio of main rotor induced downvash at tall rotor to momentum downvash 
at main rotor. A nominal value of 1.8 haa been derived from wind tunnel 
tests. Unless there are substantiating data for a particular application 
the nominal value should be used. 

589. Aircraft Gross Weight 
Expressed in lbs. 

590 . Body m 

p 

polar moment of inertia of body about body roll axis, slug-ft . ^ 

591. Body lyy 

O ' 

Polar moment of inertia of body about body pitch axis, slug- ft . 

592 . Body Izz 

2 

Polar moment of inertia of body about body yaw axis, slug-ft*. 

593. Body Ixz 

2 

Product of inertia, slug- ft . 

59 I*. Wing Option 

Fraction of O.W. carried by wing. 

595* Wing Mean Aerodynamic Chord 
Expressed in feet. 

596 . Lateral Main Rotor Shaft Tilt 
Expressed in degrees, positive tilt to left* 
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597. First Moment of Inertia Tail Rotor Blade 
Expressed in slug-feet. 

598. 

Not part of loader vector. 

599. 

Not part of loader vector. 

600. 

Not part of loader vector. 

605 - 621* Initial Values of Mode Shapes and Derivatives 

In order to reduce computer time, provision has been made to load in starting 
values of all mode shapes and their first derivatives in the event their steady 
or near steady state values are known. Values of derivatives are non-dimen- 
sionallzed by omega. 

631 - 782. 

Not part of loader vector. 

783. 

Not used 

800. Harmonics of Flatwise Bending Angles 

A +1.0 loaded into this location results in the calculation and printout of 
the harmonics of blade slope expressed in degrees at each blade station. 

801 - 836. Partial . rivatives 

In the event location U99 contains a value less than zero, the trim Iteration 
is performed using partial derivatives loaded in or previously calculated 
during prior case. Although provision has been made, for future expansion, to 
load in 36 derivatives, the present program only requires the following 9* 


Loc 813 

3Z/3A1S 

Loc 819 

3als/3AlS 

Loc 825 

9bls/dAlS 

Loc 8l>» 

3z/dBlS 

Loc 820 

3als/3BlS 

Loc 826 

obls/aBlS 

Loc 81s 


Loc 621 

3als/30^j^ 

Loc 827 




$A '.9 REV 0 



REPORT NO. SER-50910 


Sikorsky PIrcraft 


U 

OfV«lON OF iMmo 0 0<1» OPTION 

a 





It is not possi'ble to mix the modes , that is , load In some derivatives and 
have the remainder calculated internally. If the derivatives can he obtained 
by some means such as Bailey theory or derived from previous runs, this is 
highly desirable because of a substantial reduction in compute time. A word 
of caution, however. Derivatives obtained by linear theories have not proved 
reliable for high advance ratios or near stall conditions. This can result 
in divergent solutions. If the number of major iterations (Location U99) 1® 
negative and ed.1 the derivatives are loaded zero, the derivatives are auto- 
matically calculated by Bailey theory. Warning ; if successive cases are run, 
the derivatives must be reset to zero if they are to be recalculated. 

837. Gust Option 

A '*-1.0 in this location indicates that a gust penetration is to be impressed 
on the rotor after trim has been achieved. 


838. Velocity of Gust Front 

Velocity of gust front in inertial axis system. This velocity is added ko 
helicopter forward speed internally to obtain actual penetration velocity. 

See Figure 1. 

839. Gust Penetration Angle 

Angle gust front makes with rotor disc at point of tangency. Expressed in 
degrees. See Fig'ure 1. Pfige 51. 

8h0. Angle of Reference Blade 

Azimuthal position of reference blade at time of gust tangency, expressed in 
degrees. See Figure 1. 

8U1. Azimuth Print Increment in Transient 

Azimuth Increment at which stress and deflection information is to be printed 
during transient mode. Note; contents of 8Ul divided by contents of 519 Bust 
be an integer. 

8U2. Aerodynamic Print Option 

A zero loaded into this location will suppress the printing of blade aero- 
dynamic information. This only applies during the helicopter transient mode. 

8U3. Print Transient Blade Motions Option 

The printout of blade modal amplitudes and derivatives is normally suppressed 
during the transient, k *1.0 will cause a complete printout of these quanti- 
ties to occur. 
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8UU. Value of Lag Damper of Reference Blade During Transient 

In order to investigate the effect of a defective lag damper, provision for 
a different dumping coefficient of the reference blade during a tramsient 
has been provided. If this location is ignored or zero is loaded, the trial 
value is assumed. If zero is desired, a small value such a .001 should be 
loaded. 

8U5. Theta»-75 

Used in the calculation of coupled modes. Use only if different from aero- 
dynamic theta-75. 

8U6. Linear Tvist 

Used in the calculation of coupled modes. Use only if different from aero- 
dynamic linear tvist. 


Input Items Required vhen Rotor is Coupled to Flexible Support 

The following inputs are required when it is desired to couple the rotor to 
a grounded flexible support. Reference (d) describes the basis of the 
theoretical model and defines the symbols used below. 


851. Number of Flexible Support Generedized Coordinates. 5 

Generalized coordinates selected by the user to describe the support are 
typically linear and angular displacements of points on the support structure. 
The number of such coordinates is N (N £ 69). It is necessary to load H to 
define the order of the matrices describing the support structure. 

832. Not used 

853. Nximber of Modes, M 

The n\imber of support structure normal modes retained for the analysis is 
M. In general the number of normal modes must be less than or equal to the 
number of generalized coordinates (M 1 N). If the user loads non-diagonal 
mass and stiffness matrices, M will be number of support struct\ire nonnal 
modes calculated in the program, and subsequently used to reduce the support 
equations to normal form. If the user loads diagonal mass, stiffness, and 
modal matrices, M will designate the subset of the number of normal modes 
selected frem among the loaded modes that he desires to use. 
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85^ - 856. Rotor Head Linear Displacement Nxjmbers N?*, Kt* 

Rotorhead linear displacements are u, v, end w (Ref. (g)). These displace- 
ments are members of the set of generalized coordinates for the support 
structure. Each generalized coordinate is designated vith an Integer 
n\unber» these numbers ranging from one for the first generalized coordinate 
to N for the last generalized coordinate in the set. Integer numbers 
designating the generalized coordinates u, v, and v are N^*, N2*, and N3*. 
For example if there are N = 55 generalized coordinates and Nj* =20, 

N2* * 25, N3* = 51, then u, v, and w are the 20th, 25th and 51st generalized 
coordinates in the set of 55 coordinates. 

657 - 859* Rotor Head Angular Displacement Numbers, Nj, N2» N3 


The preceding remarks under 85I* - 856 apply here with N^^, N2, and N3 treated 
now as nxunbers defining the locations in the set of generalized coordinates 
of hub angular displacements 6i, 63, and 63 (Ref. (d)). 

860 . Support Str.’Cture Damping Constant, g* 

Structural damping of each normal mode of the support is represented by the 
same constant ratio g* = 2 x c/c^rit where c is the actual damping of a mode 
and ccrit is the critical damping ratio of a mode. It is emphasized that 
g* is twice the critical damping ratio c/ccRIT* 

870. 0 P(l) * 1. - Input generalized mass and stiffness matrices [m] and 

[k] 

0P(1) » 2. - Input same as for 0 P(l) except that modal matrix [yj is 
added input. 

0P(1) » 3. - Input normalized generalized mass and stiffness matrices 
[m^] and [k®3 y. See remarks in the next section referring to [y] input. 

871. 0P(2) * 0. - Input stiffness natrix [kl . 

0P(2) » 1. - Input flexibility matrix [k3” 

In the following, no printing (or action) occtirs it switches are given values 
other than those indicated. (See Reference d) 

672. 0 P( 3 ) * 1 . - Prints eigenvalues (squares of natural frequencies) and 

modal columns. 

873. 0 P(U) » 1. - Prints matrices [ m^3 » [ Cy 3 ^ [c 3 [y 3* 

87^. 0 P( 5 ) ■ 1. - Prints rotor head linear and angular accelerations and 

velocities as a function of time (rotor azimuth). 

875. 0 P( 6 ) ■ 1. - Prints modal accelerations, velocities, and amplitudea as 

a function of time (rotor azimuth). 





A 2* RCV 0 



I 


I 


Sikorsky PIrcraft 



REPORT HO. SER-50910 


876, OP(7) * 1. ~ Prints rotor head linear and angular displacements refer- 

red to undisturbed shaft axes as function of time (rotor azimuth). 

900. 0P(8) * 1. - Flexible module svitch links support module to main pro- 

gram. 

Input - Matrices 

Matrix inputs are 

[m] X N generalized mass matrix 

[k] N X N generalized stiffness ’natrix 

r Di 

L m J M X M diagonal mass matrix for normal coordinates 

Ck“] M X M diagonal mass matrix for normal coordinates 

[c3 N X N non-proportional damping matrix 

[y 3 N X M modal matrix 

The values assigned by the user to CP(l) and 0P(2) determine which of these 
matrices is input (see "Matrix Input Options"). Matrix inputs immediately 
follow loader input. F.ach matrix input is preceded by a lead card containing 
a key in colximn 6U designating the type of the matrix. Keys are 

M Mass matrix 

K Stiffness or flexibility matrix 

C Non-propoirticnal damping matrix 

G Modal matrix 
E Termination key 

A key is followied by non-zero elements of its matrix. Elements not input 
are treated as zeroes. This feature facilitates the handling of sparse matrices, 
like diagonal matrices. Null matrices require no key. Each element is pre- 
ceded by integers defining its row and colxunn subscripts respectively. Input 
of row, column, and element value is formatted 

3(13, 13, EI 5.6 

Row Col . Element 

The *3* Indicates that a aaxlmm of three matrix elements nay be placed on 
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each card. The last such non-zero element is followed by the key for the next 
matrix. This organization is repeated \aritil all non-zero elements are input. 

E in column 6U terminates these data. Only non-zero elements on and above the 
leading diagonals of [m] and [k] need to be input. The subroutine assumes 
these matrices to be symmetric and forms lover elements accordingly. How- 
ever, if the user wishes to supply the complete [ml and [k] matrices, he may 
do so vlihout penalty. When 0F(l) » 3 and [c]] is identically zero, only the 
six rows for for the six rotor head displacements need to be input. The 
elements input ere 

Yji J - Ng». N^, Ng. »3 

i - 1, 2.. M 

If OP(l) ■ 3 and [cj is not identically zero, all non-zero elements of the 
N X N matrix [y 3 are input. 

Display 

All loader and matrix card inputs are displayed to enable the user to check 
his input. Optional displays are generated according to values assigned to 
switches on input. If the switches are not assigned the values indicated 
below, printing does not occur. 

0P(3) ■ 1. - Eigenvalues (y, squares of natural frequencies) in colmn 
format, followed by each eigenvadue accompanied by its eigenvector. 

0P(U) « 1. - Diagonal stiffness and mass matrices for normal coordinates 
r k®] C m®]. Non-proportional damping in normal coordinates [y] [c] 

Lv]* 

0P(5) ■ 1. - Rotor bead linear and angular accelerations and velocities 
at each Instant of time (rotor azimuth), referred to inertial and shaft 
axes. 
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Gust front velocity in feet per second. 

V , Helicopter velocity In knots. 

H 

a , Penetration angle in degrees. 

<{/ Angle of reference blade in degrees. 

Summary of Loaded Input 

Tnere are three forms of input data; shotgun, block, and curve. The shot- 
gun loader har- the capabiU.^y of loading only values that need to be changed 
between cases. This requires the assignment of a spec'flc loader location 
for each of these items. These locations are given in numerical order In 
section three and in logical order in section four. The block loader is ^ed 
for data such at norcal modes and induced velocity vhere large blocks of data 
are loaded without the likelihood of individual numbers of a block changing 
between cases. The third type is curve data. These are in a fora that in- 
volves a dependent variable as a function of one or more independent variables. 
The prograa uses these data to Interpolate linearly between points such that 
the value of the dependent variable may be detemlned for any values tf the 
independent variables with the defined range. 
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Section 7» Common Errors and Pitfalls 

The following errors and pitfalls have been found to occur roost often in the 

operation of this and similar programs. Check carefully before running-. 

1. Improper setup of deck. 

2. Requesting incorrect options. 

3 . Leaving out word cotint on loader card. 

U. Incorrect data location. 

5 . Leaving out decimal point (when required). This will usually default 
the decimal point to the extreme right of the number field resuj.ting in 
a ridiculously large value of the parameter loaded. 

6. Delta PSI dynamic too large. This results in a mathematical instability 
of the blade dynamic response. See write up of location 518. 

7. Incorrect vsdues of the input parameters. 


Section 8. Error Messages emd Diagnostics 
Alpha Greater Than 360 

The most common error message generated by the program is "ALPHA greater 
than 360 ". This is followed by the values of the local angle of attack in 
degrees that is greater than 360 , UP, UT, PHI, THK, TWI’ST, PSI, GTQT, and 

N. 

UP, Component of velocity of air parallel to shaft, non-dimensionalized by 
tip speed. 

UT, Component of velocity of air perpendicular to shaft, non-dimensionalized 
by tip speed. 

PHI , Inflow euigle in degrees ( ARCTAN . 

THK, Blade pitch excluding twist, in radians. 

TWIST, Contribution of blade twist to angle of attack at specified station, 
in degrees. 

PSI, Azimuthal position of difficulty, expressed in radians. 

GTQT, Torsional deformation at station, expressed in radians. 
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N, Blade station of difficulty. Station one is at root. 

There is a multiplicity of causes of this error message. Examination of the 
diagonostic printout should indicate the probable cause . The following is a 
list of likely candidates. 

(a) Delta-PSI (location 518) too large. 

(b) Missing decimal point in input vector. 

(c) Improper input values. 

(d) Wrong location used in input vector. 

(e) Blade is xruly dynamically unstable. 

Word Count is MissinK 

This diagnostic can have only two causes. The word count is missing from a 
card in the data vector, or the input data are set up improperly as defined 
by the input options as requested. This could be either data out of place 
or data set up wrong. 

Program Interrupt ~ Various Information 

Although this diagnostic does not stop the calculation and can result in 
answers that may appear meaningful, under no conditions should it be ignored 
or the results believed. If the trouble cannot be found this author should 
be contacted. 

Flapping Did Not Converge 

This is a very common wexning that indicates the blade mode shapes and first 
derivatives did not converge to within the user specified tolerance after the 
prescribed maxlmian number of revolutions . Examination of the values of the 
mode shapes and their first derivatives at the 36 O degree azimuth position 
for the last and next to last revolution will reveal how serious the discrepancy 
is. The calculation will proceed after printing the warning. 

Trim Iteration Failed to Converge 

This warning indicates that either the 8 force did not converge to the heli- 
copter gross weight or the flapping did not trim to zero, within the specified 
tolerances. The final vecLues obtained axe printed for Inspection. The 
calculation proceeds after printing the warning. 
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Error in Curve Data Lookup 

In the event one of the independent parameters is out of range of the 
section data "IERR=2" is printed along with the name of the curve. The 
dependent variable is calculated using the mauc or min value of the parameter 
out of range. Calculation continues. 

In the event only one point is present on a curve the diagnostic ”IERR=5" 
is printed and the name of the ciorve. A li stout of the date should he 
checked. 

In the event a curve is requested that is not present (or the cxxrve name is 
misspelled) the name in error is printed along with all legal names. 

Section 9» General Suggestions for Running Program 

1. Do not request unnecessary mode shapes, particularly when first 
starting out with a new blade. Two flatwise, one edgewise and 
one torsional is accurate enough for initial analysis. *nie use 
of a large number of mode shapes will require a small delta PSI 
(Location 318) and result in long compute times. 

2 . Keep Delta PSI (Location 318) as large as possible, so long its 
the blade dynamics do not go unstable. 

3. When running in the transient mode, load in the partial deriva- 
tives whenever possible or use the linear derivative calculation 
option. This will reduce compute time. 

U. Use the print suppress options. This will not only reduce the 
amount of useless numbers but will reduce compute time. 

3* Always check and recheck the deck setup. Ihls is the most 

common cause of error. If the program fiM^ls, and you are unable to 
determine what has gone wrong, obtain a complete listout of the 
input before seeking help. 

Section 10. Accuracy 

The accuracy of the anedysis is best determined empirically. For any 
particular helicopter and flight condition begin with a fee modes and as 
leurge delta PSI as feasible. The number of modes can then be increased 
and the delta PSI reduced while noting the resultant change in the answers. 
When the change is smaller than the allowable accuracy the proper values 
have been established. 
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Section 11. Output Printout 
Modal Amplitudes at Beginning of Revolution 

This page of printout is designed to evaluate the convergence of the 
rotor to a steady state. It displays some data on the condition being 
analyzed, and the characteristics of the response at the beginning of eetch 
revolution . 

Output headings 6ure: 


OMEGA-R 

Rotor Tip Speed; ft/sec 

RADIUS 

Rotor Radius, ft 

RHO 

2 U 

Density of Air, lb-sec /ft 

SPEED OF SOUND 

Speed of Sound, ft/sec 

LINEAR TWIST 

Rate of Twist, the Twist 
divided by non-dimension- 
alized radius, deg. .Aerodynaia 

AIS 

LongitudineuL cyclic pitch, 
deg. 

BIS 

. Lateral cyclic pitch, deg. 

THETA-75 

Collective pitch at the 
3/ 4-radius, deg. , Aerodynamic 

LAMBDA 

Inflow ratio 

MU 

Advance ratio, V/OR, non- 
dimensional 

VELOCITY 

Aircraft forward speed, knot 


The next two lines give the flapping and lead angles at the start of the 
revoulution, deg. 


REV 


Number of revolutions of rotor 


BETA 


BETA* 


Blade root end angle for out of plane 
deflection, measured positive upward 
from precone angle, radiana 


First derivative of BETA with respect 
to azimuth angle (non-dl mens ionali zed time, 
time multiplied by rotor speed) 
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BETA** Second derivative of BETA 

DEL Blade root end angle for inplane 

deflection, measured positive for lead 
forward from prelag angle, radians 

DEL* First derivative of DEL 

DEL** Second derivative of DEL 


The next three lines give the amplitudes of the blade modes and their 
derivatives. Modes are listed from left to right, in the order in which they 
are tabulated on the previous printout page. The readings cure: 


QA Modal amplitude, non<-dimensionalized 

by radiv.8 

QA* First derivative of QA with respect to 

azimuth angle (time non-dimensionalized 
by rotor speed) , ND 

QA** Second derivative of QA with respect to 

azimuth angle, ND 


A statement at the bottom of the page will inform the user if convergence 
to a steady state is not achieved within input flapping tolerance limits. 
In any event the calculation will proceed, defining the response for the 
last revolution. 

/ 

* 

Model Amplitudes and Derivatives 

This set of output pages gives the model amplitudes and radial derivatives 
at requested azimuth angles for the last revolution. The headings are: 

PSI 

BETA 


LEAD 


Q1 to 010 

( )• 

( )•• ■ 



Azimuth angle, deg. 

Blade out of plane root end angle, positive for 
upward deflection, radians 

Blade Inplane root end angle, positive for blade 
leading, radians 

Amplitudes of blade modes 1 to 10. Zeros in 
output indicate mode is not loaded in input 

The symbols indicate first and second derivatives 
with respect to azimuth angle (non-dimensionalized 
time) of the above quantities. 
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Blade Loading and Response 


This set of output pages gives information on blade segnent angle at 
attack, airflow, euid response at each of a selected set of azimuthal locations 
for the last revolution. The data are grouped by azimuth angle which is given 
by the initial heading PSI = . The data are in columns which are tabulated 
from root to tip. The far left end right columns list the station number. 

The headings are: 


BETA 

VERT. DBF. 

lAG ANGLE 

INPLANE DEF. 
TORS. DEF. 

BR STRESS 

FLAT STRESS 

EDGE STRESS 

TORS. MOMENT 


Angle between blade segment elastic axis and plane 
perpendicular to shaft rotational axis, positive for 
deflections upward, deg. 

Out of plane deflection of center of blade segment 
referenced to blade root end, positive upward, ft. 

Angle between blade segment elastic axis and a vertical 
plane through blade root end, positive for lag aft, deg. 

Inplane deflection of blade segment, positive aft, ft. 

Elastic twist of blade segment, positive for leading 
edge up, deg. 

The sum of the flatwise bending stress and edgewise 
bending stress, psi 

Flatwise bending stress relative to blade segment 
principal axis, defined by segment bending moment and 
input value for I/C, positive for center of cxirvature 
above blade, psi. 

Edgewise bending stress relative to blade segment 
principal aucls, defined by segment bending moment and 
input value of I/C, positive for center of curvature 
ahead of blade, psi. 

Blade segment torsional moment, positive nose up at 
outboard end, in-lb 


GUST 


XCEN 


Velocity of gust at center of blade segment, par- 
allel to shaft axis, positive for cdrflov upward 
through rotor, fps or inflow velocity positive for 
upfiow fps. Depending on option selected 

Distance from center of rotation (shaft axis) to center 
of blade segnent, non-dimensionalized hy rotor radius, 
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VAR. LAMBDA 


PHI 


ALPHA 


Velocity of gust at center of blade seg^nent, par- 
allel to shaft axls» positive for airflow upward through 
rotors, fps or inflow velocity positive for upflow, f^s. 

Depending on option selected 

• • 1 / 

Tan**^ (Hp/,Ut),j the angle between the direction of air- 
flow over the blade segment and a plane perpendicular 
to the shaft, positive for upflow through the rotor 
disc and flow from leading to trailing edge of blade, 
deg. 

Aerodynrmdc angle of attack between airflow and blade 
chord line, deg. 


MACH No. 


Mach Number of airflow at blade se^ent, non-dimensional 


CL 


Lift coefficient of blade segment, non-dimensional 


CD 


Drag coefficient of blade segment, non-dimensional 


CM 


Pitching moment coefficient of blade segment, non- 
dimensional 


FW. AERO. LOAD Flatwise aerodynamic loading per spemwise distsmee 
relative to blade segment chord line, positive for 
lift on segment, Ib/in 

£W. AERO. LOAD Edgewise aerodynamic loading per spcuiwlse distance 
relative to blade segment chord line positive for 
load aft (drag) Ib/in 

TORS. AERO. )K)M. Torsional aerodynamic moment about blade segment 
quarter chord, positive nose up, in-lb/ in 


Summary of Blade Aeroelastic Behavior 

The values of various items are displayed at the indicates radial stations 
measured from the center of rotation. 

Centrifugal force in lbs 

Total max and min flatwise deflection including rigid body in feet max 
and min 

Totail max and min edgewise deflection including rigid body in feet 

Total max and min torsional deflection in degrees 

2 

Vibratory stresses in Iba/in 

2 

Average stresses in Ibs/in 
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Vibratory moments in inch-lbs (Assuming I/C's in Inch'^) 

Average moments in inch-lbs (Assuming I/C is in inch ) 

Force Integration 

This is a summary of the total rotor system loading and flight concition 
characteristics, calculated on the basis that all blades of the rotor are 
acting the ssime at a given azimuth angle. 

The first two lines are non-dimensional aerodynamic parameters based 
on aerodynamic blade loads only, defined as follows: 

2 4 

Z Force, F /np£l R 

2 2 5 

Torque Coefficient, Q/irpO R 


CZ 

CQ 

CX 

CL 

CPM 

CRM 

CPF 


2 4 

X Force, F /upQ R 
2 4 

Lift, L/irpft R 

2 5 

Pitching Moment, M /:tpfi R 
^ 2 5 

Rolling Moment, M^/irpn R 

2 4 

Propulsive Force, Fp/vpn R 
be 


These are divided by solidity o = ^ to give coefficient-solidity ratios; 

ITa 


CZ/SIGMA 

CQ/SIGMA 

CX/SIGMA 

CL/SIGMA 

CPM/SIOMA 

CPF/SIGMA 


CZ/(bc/itR) 

CQ/(bc/wR) 

CX/(bc/wR) 

CL(bc/vR) 

CPM/(bc/«R) 

CPF/(bc/i»R) 


The next two lines provide dimensional aerodynamic rotor force, moment, 
and power data. The headings are: 


Z FORCE 
X FORCE 


Rotor force F along the shaft axis, positive 
downward, lb 

Rotor force F in shaft axis system (in plane of 
rotor), positive forward, lb 




59 

PAGE 


SA 29 REV 0 








r 







I 


( 


Sikorsky PIrcraft 


U 

OtV««ON or UNfTVO C a OWPOWATOW / 

0. 


REPORT NO. SER-50910 


TORQUE 

ROLLUiG>MOM 

PITCH-MOM 

ALPHA(S) 


Y FORCE 
LIFT 

PROPULSIVE 

FORCE 

HORSEPOWER 


Rotor torque required to drive rotor, ft-lb 

Rolling moment generated by rotor in shaft axis system, 
positive for roll with right side down (right hand rule 
about forward axis) ft-lb 


Pitching moment generated by rotor, positive for pitch 
with nose up, ft-lb 


Rotor shaft angle relative to perpendicular to 


o = tan Ir 
s 


direction of flight. 

/ e: -r 

'U A 

constant inflow is used. » input value if variable 
inflow is used. 




S 


Side force F generated by rotor, positive for force 
to the right^ lb 

Lift force L generated by rotor, perpendicular to 
direction of flight, positive upward, lb 

Propulsive force F generated by rotor along 
dirrection of fli^t, positive forward, lb 

F = F sin o + F cos 

p Z 8 X ^ 

Power required to drive rotor, horsepower. 


EQUIVALENT 

DRAG 


EQUIVALENT 

L/D 


The equl'valent drag force, F^^ generated by the rotor 
along direction of flight, positive rearward, lbs 

. S?0.tHp^EPq»CTl , - 

Lift/drag ratio for the rotor, non-dimensional 
EQUIV4TEHT l/D -p- 

'ed 


EQUIVAI£NT 

PARASITE 

AREA 


The equivalent drag surface area, A^, 
rotor, as defined by 

1- V 2T 


due to the 
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The last two lines of printout provide data on the controls, flight 
condition, and sone blade properties. The headings are: 


AIS 


BIS 


Longitudinal (cosine component) cyclic pitch, positive for 
maximum angle over nose of aircraft, deg. 

Lateral (sine component) cyclic pitch, positive for maximum 
angle over left side of aircraft, deg. 


TH75(MR) Collective pitch, for main rotor, deg. 

LAMBDA Inflow ratio, non-dimensional. 

MU (MR) Advance ratio, non-dimensional 

TH75(TR) Collective pitch on tail rotor, deg. 

BODY ROLL Body roll attitude, positive right wing up. deg. 
ATTITUDE 


VELOCITY Aircraft forward velocity, knot 
BCQD/SIGMA 

A measure of the blade drag moment at a given azimuth angle, 
PSI . The magnitude shown is the leurgest value on the 
retreating half of the rotor disc. Values of O.OOU and O.OOd 
are indicative of the lower and upper stall limits (the onset 
of stall and deep stall penetration respectively). Thie 
non-dimensional parameter is defind by 


bC 


% li 

= H 

0 0 




dx 


PSI 


The blade azimuth angle at which the maximum value of 



occurs . 


BLADE WEIGHT 

Blade weight, lbs 

2 

FIRST MOM Blade first mass moment about blade root end, lb«sec 
SECOND MOMENT 

Blade second mass moment (moment of inertia) about blade root 
end, Ib-ft-sec^ 
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GAMMA 


Blade Lock Number, y> non-dimensional, defined by 
Y - PacR 


AIS FLAPPING 

Longitudinal (cosine component) of blade flapping, positive for 
blade flap up over nose, deg. 

BIS FLAPPING 

Lateral (sine con^nent) of blade flapping, positive for flap 
up over left side of aircraft, deg. 


FORCES AND MOMENTS IN BODY AXIS 

This section summarizes the forces and moments acting on the body of the 
aircraft. The individual forces and moments are given in the body axis system 
and summarized in both body and InertieQ. axis systems. The rovs of forces euid 
moments are: 

X FORCE (LBS) Positive forward 

Y FORCE (LBS) Positive toward right 

Z FORCE (lbs) Positive downward 

L MOMENT (FT-J. 1 BS) Roll moment, positive right side down 
M MOMENT (PT-LBS) Pitch moment, positive nose up 
N MOMENT (PT-LRS) Yaw moment, positive clockwise (nose right) 

The column headings 8u:e: 


ROTOR 
BODY AERO 
BODY INERTIAL 
TAIL 


Loads generated by the main rotor 
Aerodynamic loads on the body 

Inertial loading on the body, including gravity force 
Loads generated by tall rotor 


RESIDUAL Used only in the transient mode. The sum of the 

previous four elements calculated for the last steady 
state revolution before beginning the transient 
response. The residual values are subtracted from the 
four load elements in the transient condition. The 
residxial may be considered as a tolerance on steady state 
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RESICUiiL - cont'd trim. Its use prevents accelerations of the body tmtil 

a load variation occurs in the transient state. The 
residual values should be small for a true steady state 
codditioo. 


TOTALS- The toteil loaob in the body axis including residual 

BODY AXIS loads in the transient; state, excluding them in the 

steady state 


INERTIAL AXIS The total loads transformed to the inertial axes. 


SPATIAL COORDINATES 


This section describes the motions of the aircraft. For the steady state 
or trim mode of operation it displays only the forward speed CX velocity), the 
trim attitude, and sjigxilar accelerations resulting from the moment unbalances 
on the aircraft. The first line describes the following: 

PSI The azimuth angle of the first or reference blade, deg. 

GUST AT BODY Gust velocity acting on the body, positive upward, 
ft/sec. 

GUST AT TAIL Gust velocity acting on the tsdl rotor, positive 
upward, ft/sec. 


TRIM DERIVATIVES 


This section defines the changes in cyclic and collective pitch to be 
used in the next iteration in order to obtain the desired trim values for 
flapping and rotor Z force. The first matrix defines the partial derivatives 
being used as follows: 


PARTIAL MATRIX DEPENDENT ACROSS 


Aaij 

Ahi, 

AF* 

^Ig 

aaI; 

AAi, 

Aais 

Abi^ 

ABl, 

AB^* 

ABi, 

ABi, 

^»ls 

Abi^ 

AFg 

Vs 

^07 S 

Vs 


Unlt9 in degrees and pounds. 

eSi 


tr 



V 


j 
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The next oatrix is simply the inverse of the above or 


INVERSE OF PARTIALS 
ax 

(^OTE: ~ does not necessarily imply 



AAi^ 



Aai^ 

Abl^ 

AFz 




Aai^ 

Abi^ 

AFz 

^675 

\>75 

^075 

Aai^ 

obi^ 

^Z 


/; 


This matrix is multiplied by the desired changes in aig, bj^ and F^, which 
are given in that order on the next line of output headed DESIRED®CHANGES. The 
resulting changes in Aig, and 675 are given on the final line» also en- 
titled DESIRED CHANGES. For the final pass, these changes are not implemented. 
They are printed out as a measure cf the trim deviation. 

lAMOl C8 OF FLAPPING 


This section gives the harmonics of the flatwise slope cf each blade seg- 
ment, in degrees. Segment non-dimensionalized i*adius is given in the left 
column and the An and components of slope for the first ten harmonics are 
displaced. 

TRANSIENT RESPONSE 

The response of all blades in the transient mode of operation is shown. 
The azimuth angle of each blade is shown at the start with the corresponding 
response information. This is followed by a summary of loads and airframe 
response. The output parameters are all cs previously described. 

TRANSIENT RESPONSE SUMMARY 


The detailed transient response is followed by a sussnary showing the varia- 
tion of key parameters with time, when time is defined by the azimuth angle of 
reference hlade. The columns show for the reference blade: 

1) Azimuth angle, deg. 

2) Root flatwise stress, psi 
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TRANSIENT RESPONSE SUMMARY (cont'd) 

3) Flatwise stress at station 8, psi 
U) Root torsion moment, in-lb 

5) Rotor Z force, lb 

6) Body aerodynamic Z force, lb 

7) Total Z forci on the body, lb 

8) X force on the aircraft, inertial axes, lb 

9) Y force on the aircraft, inertial ax^ss, lb 

10) Z force on the aircraft, inertial axes, lb 

11) Roll moment on the aircraft, Inertial axes, ft-lb 

12) Pitch moment or. the aircraft, inertial axes, ft-lb 

13) Yaw moment on the aircraft, inertial axes, ft-lb 
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i.ne folloving page contains the definitions ,f the important output quantities 
and their associated sign conventions. 
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onwON V 



OCM^OiiATION 


APPENDIX A 


A listcut of a complete run including the "JCL” for each of the 
three steps follows. While the specific "JCL’* may vary from installation 
to installation the basic flow will remain the same. The particiilar **JCL” 
presented here is set up to run each of the steps automatically as one 
batch submission. An alternate fonn could be constructed to run each step 
iniividualiy . The former has the advantage of less elapsed time, the 
Ivitter rllovs for scrutiny of each individual step before proceeding to the 
next^ 
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Additional Input for Unsteady Aerodynamics 


The unsteady curve option requires the following peripheral curves in addition 
to the basic unsteady data. 


Curve 


ASTCN 

RCNM 

ASTCM 

ClIREGl 

CNREG2 

CMREGl 

CI'1REG2 


Cxirve 

Title 


Value of stall alpha associated with CN as function of Mach No. 

Ratio of CN stsdl at actual Mach no. to CN stall at Mach no. = 0. 

Value of stall alpha associated with CM as function of Mach no. 

Upper boundary of alpha associated with EN as a function of 

Mach no. to define end of unsteady region. 

liower boundary of alpha associated with EN as a function of 
Mach no. to define beginning of steady region. 

Same comment as CNREGl except for CM. 

Same comment as CNREG 2 except for CM. 


The area between the upper boundary of the unsteady data and the lower 
boundary of the steady data is a transition region. If the o, Itech number 
conbinatior. is in the region a linear interpolation between unsteady and 
steady is performed to calculate CN or CM. 

• - 

In addition the maximum and rainimiun values of A and B( functions of ® ®nd a) 
are defined in a data statement in subroutine "AEROD”. Thus if these values 
are to be changed the subroutine "AEROD" must be recempiled. A complete 
description of the mechanics of the unsteady aerodynamics may be found in 
reference (e). 


AdditionaJ. Input for Induced Velocity Calculation 

The following three cui'ves are used by the aeroelastlc module to prepare 
data files for the calculation of Irduced velocities. See Appendix C. for 
a complete description. 

Curve Curve 

Title 

ACLZ Alpha for zero Cj^ vs. Mach No. 

ACLM Alpha for max vs. Mach No. 

ASLP Lift curve slope vs. Mach No. 
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A sample listout of all necessary data assuming a 0012 airfoil 
follows. The "JCL" to set this data up on File 22 is also included. 
However this data may he loaded directly into the program via the card 
reader if desired; see location 135 on page 19. 1 
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TECHNICAL APPROACH 
General Description and Assumptions 

The UARL Prescribed Rotor Wake Inflow Program (hereafter referred to as 
the rotor inflow program) computes the circulation and inflow distributions 
along the blades based a prescribed wake model and a number of assumptions 
regarding tht ierodyn:- . cs at each blade section. The section operating 
conditions a.re t-^enerally prescribed from blade motion and control information 
obtained fror a separate blade response program. The wake model is either 
generated internally in the con^puter program (undistorted wake) or prescribed 
from a separate program. Procedurrr; have been developed for coupling the 
programs together in a linked manner as shown in Figure 1. Descriptions and 
results pertaining to the development and application of the rotor inflow 
program in combination with wake geometry and blade response programs are 
presented in References 1 through U. 

Input to the rotor inflow program is obtained from the coupled mode 
dynamic analysis. The latter program is often referred to as deck or program 
Yl4l in this appendix. This YlUl should not be confused with the uncoupled 
aeroelastic program described in Reference d of the text. The YlUl referred 
to in this appendix is the coupled mode dynamic analysis described in 
Reference c. 

Briefly, the mathematical model in the rotor inflow program consists of 
the representation of each blade by a segmented lifting line, and the helical 
wake of the rotor by discrete segmented vortex filaments consisting of 
trailing vorticity which result from the spanwise variation of bound circula- 
tion. The circulation of the wake for each blade changes with azimuth position 
auid is periodic for each rotor re\rolution. The blades are divided into a 
finite number of radial segments, and the induced velocity at the center of 
t »ch selected blade segment is computed by summing the contributions of each 
bound and trailing wake segment. 'The contribution of each vortex segment is 
obtained through use of the Biot-Savart equation which expresses the induced 
velocity in terms of the circulation strength of the vortex segment euid its 
geometrical position relative to the blade segment at which the induced 
velocity is desired. The bound circalation distribution is determined by 
relating the wake circtilations to the bound circulations, expressing the wake 
induced velocities in terms of the unknown bound vortex strengths by means of 
the Biot-Savart equation, and developing a set of simultaneous equations 
relating the bound circulation and local blade angle of attack at each blade 
segment. These equations thus involve the known flight condition, wake 
geometry, lift-curve slope, and blade motion and control parameters and the 
unknown botind circulation values. Solution of these equations yields the 
desired bound circulation values, which, when combined with the appropriate 
geometrical relations in the Biot-Savart law, produce the required induced 
velocity distribution. 


^0 

I 

V 



PAGE 



SER-50910 



4^ O 
O t-i u 
fXi (Q fU 




SER-50910 


Figure 2 



Sequence of Frogrea Operations 
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A simplified flov diagram shoving the sequence of operations of the 
program Is presented In Figure 2, The section parameters (airfoil character^* 
istlcs and geometric angle of attack) are computed. The hlade-wake geometry 
is formulated. The wake Influence coefficients (geometric coefficient) at 
the blades are calculated using the Biot-Savart law. Numerical techniques 
6u*e used to solve the circulation matrix. Finally, the inflow (induced 
velocities) Eure comrated. 

The fundamental background for the analyticEil fonatilation used in the 
computer program was obtained from the technical approach described by 
Plziali and DuWaldt in Reference 3- Although a new computer program has been 
developed, and many modifications and refinements have been incorporated, 
many of the fundamental assumptions are retained from that reference. 

The rotor inflow program is based on the following formulation and 
assumptions: 

1. Each blade is represented by a lifting line (bound vortex) divided 
into a finite number of segments (blade segments) each having a different 
circulation strength (see Figure 3). The aerodynamic chenracteristics at 
the centers of each segment are ass^uned to be representative of the entire 
segment over a finite azimuth interval. The use of lifting line rather than 
lifting surface theory results in the satisfying of the aerodynamic boundary 
conditions at only one point on the chord (qxiarter chord point) of each blade 
segment . 






2. The wake is represented by a finite nxsnber of vortex filaments 
trailing from the blade segment boundaries. Each filament is divided Into 
straight segments, the lengths of which are determined by a specified wake 
azimuth interval which is equivalent to the azimuth Interval of each blade 
(see Figure 3). The circulation strength of each trailing vortex segment is 
constant along its length, and is equivadent to the difference between the 
circulation values of its culjacent bound vortex segments in accordance with 
the Helmholtz laws of conservation of vorticity. The circulation strengths 
of different vortex segments along a vortex filament vary in accordance 
with the variation of the bound circulations with azimuth position. 

3. Viscous dissipation effects on the wake circulation strengths are 
neglected in that the circulation of a given wake segment is constant with 
time. However, the number of wake revolutions retained in the analysis can 
be limited to evaluate an abrupt dissipation of the wake. Also, the viscous 
rollup of the tip vortex can be approximated by combining tip filaments 
beyond a prescribed rollup aximuth Interval. In addition, a vortex core size 
Is assigned to each vortex filament. Potential theory Is assumed to apply 
outside the vortex core. Inside the vortex core zero flow Is assumed. 

k» It Is assumed that the rotor Is operating In steady-state fll^t. 
The Inflow and wake from each blade is assumed to be periodic with blade 
spacing. This Is, the inflow and wake geometry Is the same for each blade 
when at a specific azimuth position. 
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Figure 3. Representation of Blade and Wake by 
Bound and Trailing Vortex Segments 
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5. Shed wake segments (segments normad to the trailing vake segments 
mentioned above) arising from the time rather than radied variations of 
boxind vorticity, are not included in the watke model. Although this omission 
technically violates the Helmholtz lav, it is believed that a more acciirate 
representation of the shed wake effects is obtained through the use of 
experimental unsteady airfoil characteristics in the analysis. This implies 
that the primary effects of the shed vorticity aure those associated with the 
shed vaike immediately behind the blade, and thus cam be approximated by those 
of a fixed-wing type of wake. Miller (Reference 6 ) shows that this is 
reasonable at rotor advance ratios usually of Interest. This approach not only 
permits a factor of two reduction in computer time amd prevents unreadistic 
results associated with a finite filament wake model, but more importamtly, 
permits nonlinear unsteady stall effects to be included in a rational manner 

in the blade response program. 

6 . The wake geometry is prescribed from amalyticad or experimentad 
res\dts. Various options are available for selecting the representation of 
the wake model. 

7. The airfoil at the blade is assumed two-dimensional (radial velocity 
components are neglected). For the lineaurized circulation solutions and 
associated cedculations , a set of lift curve slopes, stall angles of attack, 
amd angles of zero lift are provided which vary with Mach number. These may 
be based on unsteady aerodynamic data, if available, and provided directly 
from a blade response program. Below stall the lift curve slope is amsumed 
constant. The blade section circulation is limited to a constamt value for 
angles of attack above stall. 

8 . In the blade-wake geom'^try calculations the blaules 80*0 amsumed 
straight (rigid blades). However, flexibility effects may be included in the 
circulation solution by providing the necessaury noninduced velocity at eamh 
segment associated with flexible blade motions and controls from a blade 
response program. Also, an anhedral tip may be prescribed. 

9. Tangential Induced velocity components are neglected. 

10 . angle auisumptiona are included in the circulation solution. 

11 . The aerodynamic interference effects of the rotor hub, fuselage, 

etc. are neglected. The rotor is assumed to be operating out-of-ground effect. 
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Wake Geometry 


The rotor inflow program requires that the rotor wake geometry be speci- 
fied in order for circulations and induced velocities to be determined. There 
are several alternatives for rotor wake geometry. The least complex wake 
geometry is a classical undistorted wake which is simply a function of the 
flight condition and mementum inflow velocity. The coordinates of this 
helical wake, which is skewed in forware flight, are easily generated in a 
prescribed wake type of analysis given the governing parameters which can be 
iterated on if desired. If analytical, realistic distorted wake geometries 
require more complex and operationally expensive computer analyses. Appro- 
priate experimental wake data is certainly most desirable, but is not currently 
available for most rotor configurations and fliglit conditions. It has been 
established that the requirement for distortions from the classical type 
wake geometry is dependent on the rotorcraft configuration, flij^t condition, 
parameter of interest and accuracy required. For example, an iindistorted wake 
geometry is generally sufficiently accurate for integrated performance calcula- 
tions (thrust, torque, etc.) for conventional rotorcraft operating in hi^ 
speed flight. However, low speed flight conditions wake distortion? are very 
significant . 

Complete generality regarding the specification of the wake geometry is 
provided in the computer program. This was accomplished by requiring that the 
wake coordinates for the wake segment end points be stored on tape for compu- 
tational purposes in a prescribed format, thus no requirement was made as to 
how these coordinates are obtained. This facilitates the adaption of improved 
wake models to the program as they become available. For this program at the 
present time several options are available for wake geometry. The following 
three types of wake models are presently available: (l) A classical undistorted 

wake model option is available internally, which computes wake coordinates and 
stores them on tape, (2) an analytical weke model (undistorted or distorted 
wake) can be obtained from a separate wake geometry program which can be run 
linked to this program, and (3) a provision is available for a classical wake 
with a tip filament overwrite using available experimental data. 
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Undlstorted Wake Model - In its sin^jlest form, the wake can he assumed 
to be a classical uadistorted skeved helical sheet of vorticity defined 
from momentiam considerations (hereafter referred to as the laidistorted wake). 

A sample undistorted wake representation is shown in Fig. 1*. The coordinates 
of the undistorted wake representation sure obtained from the rotor advemce 
ratio, , tlirust coefficient. Op, and smgle of attack relative to the tip 
path plane, tPP* example, the top view of the tip vertex filaments 

is obtained directly from the helicoidal path of the bla<1e tip as it translates 
at the velocity V^os TPP rotates at the velocity R. The side view is 
dependent on the wake skew angle (angle between normeu to the rotor disc 
and wake boundary), which is normally defined from momentum considerations. 

Subroutine CLASWK is used to provide a classical undistorted wake 
geometry. The axial transport velocity of all wake segment end i>oints Is 
equal to a constant which is an input item. For example, the momentum 
induced velocity has norm6U.ly been used for this value. The surface boundary 
of the resulting wake model is therefore a skewed cylinder in shape. Each 
wake filament thus forms a skewed helix. 


Distorted Wake Model - A distorted wake model may be incorporated 
in the analysis via tape input. A rotor wake geometiy analysis, such as 
that developed at UARL and described in References 1 and 2, can be iised 
to generate wake coordinates on tape for use in the rotor inflow analysis. 

In addition to tape input, a provision is available in the inflow 
program for overwriting the coordinates of the undistorted tip filaments, 
generated internally, with tip vortex coordinates determined from another 
source (ejqperimental or analytical) . This option can be used to more 
accurately model the important tip vortex geometry. 
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Wake Coordinate System - A Cartesian coordinate system is used in the 
program to define the wake geometry. The coordinate system is defined as 
follows : 

1) Coordinates arc nondimensionalized by the rotor rculius. 

2 ) Origin lies at the center of the rotor hub. 

3) The x-y plane is parallel to the rotor tip path plane. 

U) The coordinate system is a right handed system with z positive upwards 

5) The coordinate system is aligned such that positive x is always in 

the downstream direction and zero azimuth lies on the positive x 
axis. 


Circulation and Induced Velocity Solution 

Following the formulation of the wake geometry representation in the 
program, the geometric relations between the wake and the blades are calculated 
and the equations for the blade circulations and induced velocities are applied 
The Kutta-Joukowski law and the Blot-Savart Inw are the basic relations used 
to obtain a closed form solution for these parameters. The Kutta-Joxikowskl 
law relates the blade circulation and induced velocity distributions. The 
BiotoSavart law relates the blade induced velocity distribution to the wake 
circulations and wake geometry. The term "geometric coefficient" is used to 
describe the influence coefficients in the latter relation, which are functions 
only of the wake geometry. Several assumptions, mentioned above, for the 
airfoil lift characteristics, velocity components, and inflow euigles were 
Included in the establishment of the circulation matrix in order to linearize 
the solution. The solution procedure consists mainly of (l) determination of 
the geometric coefficients, (2) establishment of the circulation matrix, 

(3) solution of the circulation matrix by a modified Gauss-Seidel iteration 
technique, (U) solution of the induced velocity distribution (inflow) using 
the results of the circulation solution, and (5) harmonic analysis of the 
inflow solution and calculation of associated parameters. 

In order to summarize the technical approach and Indicate where 
assumptions are made, the primar>' equations contained in the circulation and 
inflow solution and a brief derivation thereof are presented below. 
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In order to relate to the program input-output symbols and to describe 
some of the symbols generated within the program, program symbol notation is 
used in the following equations. 


The blade circulation at any blade segment I, for any azimuthal position 
J, is; from the Kutta-Joukovski relation. 


where 


CIRC(I,J) » C(I) • U(I,J) • CL(I,J)/2 


CIRC(I,J) 

C(I) 

U(I,J) 

CL(I.J) 


2 

^ local blade element circulation, ft /sec 
= local blade element chord, ft 
» local blade element resultant velocity, fps 
* local blade element lift coefficient 


( 1 ) 


The local blade element resultant velocity U(I,J) can be expressed as 

U(I,J)^ « [UT(I,J) ♦ UTI(I,J)] ^ + [wi(l,j) ♦ W(I,J)] ^ (2) 

where UT(I,J) » OMGR . (RS(I) + MU*Sin(PSI)) (3) 

and OMGR = Rotor tip speed, fps 

MU » Rotor advance ratio 

RS(I) « Rotor blade segment radial position (nondimensionslized 
by R) 

UTI(I,J) ■ Induced tangential velocity, fps 
W(I,J) « Induced inflow velocity, fps 


and with small angle assumptions 

Ui(I,J) ■ Noninduced inflow velocity (flapping, pitching, aeroelastics , 
etc.), fps. 

For helicopter rotors, UTI(I,J) may be neglected. Also since W(X,J) and Wl(l,J) 
are small in the important blade tip region relative to the rotational velocity, 
it is assumed that 


U(I,J) - UT(I,J) 


(U) 
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SO that a linear closed form solution can be obtained. Also, in order to 
linearize the solution, the lift coefficient is expressed as 

CL(I,J) » AT(T,J) • ALPHAE(I,J) (5) 

lAere AT(I,J) » local blade element lift curve slope 

ALPHAE(I,J) » local blade element angle of at*’*ok 

In stall, the maximum lift coefficient is limited in the program to CLMAX(I,w; 
which is determined from input values of angle of attack values (ALMAX(I,J) 
beyond which stall is assumed. The angle of attack is, for small angle 
assumptions and consistent with the above assumptions, 

ALPHAE(I,J) » [wi(l,J) + W(I,J)] /UT(I,J) (6) 

Substituting Equations {k) through (6) into Equation (l), the circulation 
equation becomes, 

CIRC(1,J) » C(I) • AT(I,J) • [wi(I,J) + W(I,J)] /2 (7) 

Since there are MSIZE unknowns, 

where M5I7.B * (ITOT +1) • JTOT and ITOT and JTOT are the total number 

of blade segments and azimuth positions, 
there are ^BIZE eonations of the form of Equation (7). This set of fCIZE 
equations can be written in matrix notation as: 

C1RC(I,J) - [c(D • AT(I,J) • {wi(I,J) ♦ W(1,J)| ] (8) 

The axial induced velocity in Equations (7) and (8) is expressed, using 
the Biot-Savart law, as a function of the trailing wake filaments and circula- 
tion strengths. The axial velocity Induced by a given wake segment (subscript 
M) of filament K at a given blade segment (I,J) is expressed in teraw of the 
Blot-Savart law as 

W(K,M) ■ CIRCT(K,M)/l|pR • GC(K,M) (9) 

idiere CIRCT(K,M) • circulation strength of the Hth segment of the Kth 

filament, ft^/sec 

GC(K,M) " geometric influence coefficient from the Biot-Savart lav 
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The geometric coefficient, GC, is a function only of the prescribed wake 
coordinates (X(K,M), Y(K,M), Z(K,M) . The lengthy expressions for CC may be 
obtained from the program listing included herein (Subroutine Jack). Since 
the circulation distribution is assumed to be periodic in time for the forward 
flight condition, each wake filament has a periodic circulation strength 
distribution along its length. The total axial induced velocity due to 
contributions of all trailing circulations for all of the filaments at the 
blade point (I.J) is the sum of Equation ( 9 ) evaluated for all segments in the 
wake. The resultant expression from Equation (9) f comes 

KL.IM MLIM 

W(I.J) » I t CIBCT(K.M) • GC(K,M) (lO) 

MirR g M 


where KLIM is the total number of trailing filaments, and MLTM is the total 
number of segment end points for esch filament. 


However since the circulations are periodic in M this relationship can be 
expressed by combining elements with like clrcxilations, 

KLIM JTOT 

Z V CIBCT(K,MM) * GCTBOW(K,MH) 

. ^ 

K MM (11) 


ehere 

GCTBOW(K,MM) 


MLIM 

(M>MM) 


GC(K,M) 


MM combines all elements of a vortex filament shed from a blade at 
the same blade azimuth position, and thus combines elements with like 
circulations. The total induced velocity due to all bound circulations of 
the rotor can be expressed as, 

B ITOT JTOT 

WB(1,J) - 1 Z I Z CIRC(X,MM) * GCB(I,MM) (12) 

kirR IB 1 


where CCB(I,MM) • geometric coefficient due to bound circulation (using 

Biot-Savart lav) 

VB(I,J) ■ induced axial velocity due to bound circulpLion 

IB ■ nui^er of blades 

ITOT ■ number of bound segments per blade 


1 
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In order to express the total axial induced velocity due to all of the trailing 
circulations in terms of the bound circulation, the circulation strength of 
each trailing filament segment may be expressed as the difference between the 
bound circulation of the two blade segments adjacent to its point of origin 
at the blade at the time of generation. Thus, 

CIRCT(K,MM) « CIRC(K,MM) - CIRC(K-1,MM) (13> 


Substituting E(iuation (12) and Equation (13) in Equation (ll), the total 
axial induced velocity becomes 

1 KLIM JTOT 

y(l^j) + :: — J Y. CIRCT(K,MM) • GCTROtf(K,MM) 

1»#R K MM 

1 B ITOT JTOT 

+ Z Z Z CIRC(I,M) • GCB(I,MM) (lU) 

l»rR IB I MM 

Rearranging Equation (l**) to combine like circulations leads to 


where 


1 B ITOT 

w(i,j) - — z z 

it«R IB II 

GCDIM(II,MM) * GCTR0W(II+1 


JTOT 

1 

MM 


,MM) - 


CIRC(II,MM) • 
GCTROW(II,MM) 


GCDIM(II,MM) 


(15) 


+ GCT(II,MM) 


Sxibstitution of Equation (15) in Equation ( 8 ) results in a mtrix equation 
in which the only \ujknovns are the bound circulation, CIRCU»Ji 

[giro] « [const] ♦ ^gcmat} [giro] 

where GCMAT > final matrix influence coefficients 

CONST ■ matrix constants which are only functions of the known 
input quantities. 

Solving for CIRC: . 

[ciRc] ■ [cONST^ |i-OCMAt} UT) 

A Gauss*^eidal numerical iteration technique for evaluating this circulation 
matrix is included in the computer program. Following the circulation solution, 
the induced velocity distribution is obtained from Equation (8). 
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COMPUTER SYSTEM INFORMATION AND PROC®AM STRUCTURE 


The UARIj Prescribed Rotor Wake Inflow Analysis Program (Deck P 389 ) Is 
programmed in Fortran V computer Ismguage for use on the UNIVAC 1110 hi^> 
speed digital computer. The program has been converted for use on IBM 370 and 
360 computers. 

The structure of the program consists of an initial routine (F3d9)» a 
major subroutine (JACK), and a series of minor subroutines. The symbolic 
names for the progi'sun subroutines are listed below. A brief description 
of each of the subroutines is presented in Appendix I. A diagram of the 
sub. outine calling sequence is presented in Figure 5. 


AERPAR 

ANALYS 

CIRPRT 

CLASWK 

FLAPIT 

F 389 

IRTERP 

INTPO 

JACK 

LINKUP 

MLS0P2 

PCOORD 

PRINT2 

REW 

TIPWAK 

TTRAP 


PR0GRA.M SUBROUTINES 


BLAGE 

BLDXYZ 

DATAIN 

DISVAK 

HARM 

INDAT 

INTPOE 

IHTVEL 

LOADER 

MILSOP 

PLOTrr 

PRINIT 

RREC 

SKPnL 

UAEBO 

WE7 


WREC 
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Figure 5. Subioutine Calling Sequence Diagrav 
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SUBROUTINE CIASWK 

This subroutine computes a classical undistorted vake. The vake 
is hub referenced and parallel to the tip path plane. After the wake 
has been generated it is written on tape. If blade lagging is desired. 
Subroutine BLAGE is called. Subroutine £EW, W£F and WREC are also called 
from this subroutine. 


SUBROUTINE DATAIN 

This subroutine controls input of basic variable data and computes 
control associated parameters and constar.ts. These parameters and constants 
are stored in named common blocks to be used throughout the whole program. 

It calls subroutine LOADER. 


SUBROUTINE DISWAK 

This subroutine controls the input wake geometry from tape cmd the 
associated wake print option. It is called for each rotor position. 

The wake geometry is read from tape by a tape reading subroutine SREC. 
The tip wake overwrite option is used in this subroutine. The following 
subroutines are called from this program: PCOORD, REW, SKPFIL, TIFWAK. 


SUBROUTINE FLAPIT 

This subroutine computes or reads in by cards the blade element 
noninduced axial velocity component (Wl). Rigid blade equations are used if 
W1 is internally computed. Wl can be input by cards via this subroutine. 


SUBROUTINE F38o 

This is the m-*n program and controls the flow of the program. The 
following subroutines are called in F 389 : ANALYS, BLDXYZ, CIRPRT, CLASWK, 
DATAIN, INTVEL ^ACK, MILSOP, MILS0P2, PRINIT, TTRAP. 
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SUBROUTINE HARM 

This subroutine hannonically analyzes data and computes the hannoaie 
coefficients for a positive Fourier series. 


SUBROUTINE INDAT 

This subroutine is used when running linked vith another program. 
It loads into the variable input array a set of standsurd input items. 


SUBROUTINE INTER? 

This subroutine is a linear interpolation subroutine. It is used 
throughout the program wherever linear interpolation is required. 

SUBROUTINE INTPO 

This is a print option subroutine which used fixed point output 
format. It is used for printing out variables sis a function of blade 
aximuth and blade radial position. 


SUBROUTINE IHTPOE 

This is a print option subroutine which uses an exponential output 
format. It is also used to output information as a function of blade 
azimuth position and blade radial station. 


SUBROUTINE INTVEL 

This subroutine is used to control the interpolation of infloif 
velocities and the output of these velocities. This subroutine calls 
subroutine INTER?. 


SUBROUTINZ JACK 

This subroutine is the major ccxnputational subroutine. It computes 
the geometric influence coefficients of each wake element at each blade 
r ad ial station and azimuth position. Geometric coefficients for both bound 
and trailing vortex elements are computed in this subroutine. It calls 
subroutines A£R?AR, DISWAK. FLAPIT. ?RIRT2, UAERO, UREC. 
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SUBROUTINE LINKUP 

This subroutine links another program to this program and allovs 
transfered of data via drum. 


SUBROUTINE LOADER 

This subroutine reads variable input data from cards with a 
specified format. 


SUBROUTINE MILSOP 

This subroutine solves for circulations using a Gauss Seidel iteration 
technique. If the circulation matrix Is "small" (l.e., can be stored in 
the computer), the geometric influence coefficients are read from drum 
into core and stored In matrix form. The above technique is then applied 
to determine the circulation solutions. If the circulation matrix is 
large Subroutine I€jS0P2 Is used. Incorporated Into the iteration technique 
Is a check on maximum circulation as determined from stall considerations. 

If cuiy guess In the Iteration exceeds the maximum circulation, this guess 
Is reset to the maximum clrcxilatlon and iteration continues. It calls 
Subroutines REN, and RREC. 


SUBROUTINE MILS0P2 

This subroutine solves for circulations using tne Gauss Seidel 
technique when the matrix Is "large" (l.e., matrix is too large to store 
In the computer). The only difference between this subroutine and the 
previous subroutine MILSOP is that the matrix Is manipulated on drum during 
the iteration sequence (matrix is not stored in core). It calls subroutines 
REM, RREC, «nd MRBC. 


SUBROUTINE PCOORD 

This subroutine Is used to print wake geometry coordinates. 


SUBROUTINE PRINIT 

This subroutine prints the basic variable Input data In a special 
output format which is compatible with the loader Input requirements. 

It also prints the variable's location in the arrsy. name description, 
units and value. ^ \ 
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SUBROUriME REW 

This subroutine rewinds a desired tape or drum unit when called. 
SUBROUTINE RREC 

% 

This subroutine reads a record of specified length from drum or tape. 
SUBROUTINE SKPFIL 

This subroutine skips a file on tape when called. 

SUBROUTINE TIFWAK 

This subroutine will read tip wake filament geometry from cards 
when called, and overwrites the classical wake tip filament geometry- 


SUBROUTINE TTRAP' 

This is a time trap subroutine idiich confutes the time required to 
complete various portions of the program. 


SUBROUTINE UAERO 

This subroutine confutes sectional airfoil characteristics from 
information obtained in card form from another program. It uses 
subroutine INTER? . 


SUBROUTINE WEF 

This 8'd>routine will write an "end of file" on either drum or 
tape iriien called. 


SUBROUTINE WREC 

This subroutine is xised for drum writing purposes. It writes • 
record of specified length on either drum or tape. 
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Modes of Operation (Linked vs. Independent) 


The rotor inflow program (Deck F389) can he used In either of two modes 
of operation •— linked or independent. 

For the linked mode of operation, the inflow program is coupled (linked) 
directly to the hlade response program, Yll*l (see Fig. 1) and the two 
programs are nui as one. That is the blade response program provides the 
necessary input to the inflow program, and the inflow program then computes 
and transmits the Inflow solutions to the blade response program. In this 
manner of operation, unless overwritten with separate input for Deck F389* 
all of the input items necessary for execution of the Inflow program are 
transmitted from the blade response program or preset internally in the 
Inflow program. The flight condition, rotor design, rotor pitch control, 
blade response, and airfoil data parameters sure input to or computed in 
the blade response program (Yl^l) and transmitted internally to the inflow 
program (F369}. Most of the rotor and wake model parameters and option 
controls are preset internally in the inflow program. This relieves the 
requirement for input to the inflow program for linked operation other than 
two control cards to be described. However, caution should be exercised 
by reviewing the input instructions to Insure that the transmitted and 
preset values conform with the requirements for the specific rotor configure* 
tion and operating condition of interest. Any of the transmitted or preset 
values may be overwritten, if desired, with the appropriate selection of 
input items via a secondary LOADER input provision to be discussed. 

Generally, it is necessary to overwrite at least a few items. For some 
cases selection of certain items (e.g., inflow stations, wake roU*up, core 
radius, etc.) may be critical for accurate solutions. Thus it is 
recommended that the user become familiar with all of the Infl'W Itrogram 
input items for effective use of the program. 

The second mode of operation of the Inflow program is the independent 
(stand alone) mode, whereby the program functions without being coupled 
to the blade response program. For this mode of operation, all uf the 
input items must be input via cards. 
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DESCRIPTION OF PROGRAM INPUT 
SijmEiary of - Input Parameters 


The input required for the rotor inflow program are summarized below 
by grouping the input parameters into descriptive categories. 

Flight Condition Parameters; Flight velocity, rotor tip speed, 

rotor shaft angle, speed of sound 

t- 

Rotor Design Parameters: Number of blades, flapping hinge, 

offset ratio, rotor radius, blade 
chord, blade twist, blade precone 
(for hingeless rotor), tip anhedral 

Rotor Pitch Control Pararieters: Collective pitch, cyclic pitch 

Blade Response Parameters: Coning (articulated rotor), flapping 

angles, lag angle, blade bending and 
torsionad flexibility (optional) 

Rotor Model Parameters: Radiad stations, azimuth increment 

Wake Model Paurameters: Azimuth increment, number of vaOce 

revolutions, wake transport velocity, 
vortex core radius, number of vortex 
filaments to rollup in tip vortex, tip 
vortex rollup azimuth, inboard wake 
truncation azimuth, tip vortex 
coordinates (optionad) 

Airfoil Data: Lift curve slopes, stall amgles of ■ 

attack, angles of attack for zero lift 
(all vs. Mach number) 

Option Controls: Print, punch, link, auad input controls 

Although the r-vnber of input parameters appears large, provisions 
have been incoirporated in the program to facilitate the input setup by: 

( 1 ) omitting input for unnecessary parameters for specific cases, 

(2) internally prescribing input quantities when operating in the 

I linked mode, and 

(3) avoiding input repetition fbr mult.lpl« cases. 
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Sequence of Input 


The input sequence for the computer ctird input is: 

1. Basic Vcuriable Input Data (LOADER FORMAT) 

2* Multiple Case Control (1 Card) 

3. Secondary Changes in Variable Data, LOADER FORMAT, (Optional) 

4. Multiple Case Control (1 Card) 

5 . Sectional Airfoil Characteristics Input (Optional) 

6. W1 Input (Optional) 

7. Tip Vortex Coordinate Overwrite Input (Optional) 

For runs with isultiple eases, the input sequence for eases folloving 
the first case is identical except that repetition of data in the Basie 
Variable Data is not required.* , 

A detailed description of the input items follows below. A sasple 
listing of the input data in computer card format and computer printout 
format is presented in a later section of this report. 


* A eoi?puter ease is referred to as the inflow solution for one 

CMDbination of rotor design, operating condition, and wake geometry. 
A computer run is referred to as the inflow solution for one ease or 
multiple eases submitted to the computer at one time. 
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Deck Deck in 

F389 YlUl* Deck F389 DESCRIPTION OF INPUT ITEM 


1 517 V 

2 1 OMGR 

3 l»66 SOUND 


Flight velocity, kt. This item is used in 
computations involving various blade element 
aerodynamic quantities. 

Rotor tip speed, fps. Product of the rotor 
rotational velocity (rad/sec) and the rotor 
radius (ft). 

Speed of sound, fps. This item is used in 
calculation of local blade element Mach 
number . 


U DPSI Azimuth increment, deg. This item defines 

the computational azimuth increment of a blade 
over one revolution. It also is used to 
divide the trailing vake filaments into 
elements and to define the azimuth interval 
of the circulation and inflow solution. 

360/ (B* DPSI} must be an integer. DPSI 
must be equal to or greater than 10 deg. 15 
deg is normally used for ?, 3, It, 6, or 
8 blades, although 30 deg is used when 
computer cost is more of a concern than 
increased accuracy. (Linked Version: 
prescribed as 15 deg except for five blades 
(12 deg) and seven blades (17.1^3 deg). 
Transferred from Deck YlUl). 


5 2 R 


Rotor radius, ft. 


6 130 B Number of blades 

7_21 C Blade chord for each segment in the inflow 

solution, ft. Maximum of 15. Ordered from 
the innermost of the outermost segment. 
(Linked Version: interpolated internally at 

radial stations for inflow solution(locations 
2U-38) from input chord lengths in locations 
150 >i 6U of Deck YlUl). 

• LiiAed Version. Locations in Deck Yll»l for automatic transfer of Input 
to appropriate locations in F3d9 module. 
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DESCRIPTION OF IRPOT ITEM 


22 


23 


E 


STN8 


2U.38 


RS 


Blade flapping hinge offset ratio. Offset 
nondimensionalized by the rotor radius* R. 

Used to obtain blade segment orientation. 

Number of blade segments used for the Inflow 
solution. A maximum of 13 is cdlowed. 

Normally, 9 se^sents are used. (Linked 
Version: prescribed internally as 9.) 

Radial coordinates of blade segment centers 
(blade stations) for the inflow solution 
nondimensionalized by the rotor radius, R. 

These coordinates are used to define the 
representative locations of the blade 
se^nents for which the aerodynamic parameters 
are calculated and the inflow solutions are 
obtained. In addition, the boundaries of the 
blade segments, which designate the origins 
of the trailing vortex filaments, are 
calculated from these segment center locations. 
The Innermost sepaent center is normally 
selected so as to position the corresponding 
innermost segaent boundary at the blade root 
cutout position (inner aerodynamic boundary 
of blade). The maximta number of blade 
se0Dents is 13. Careful selection of the R8 
values must be sade to insure that the radial 
position of each segment boundary between 
adjacent segments is Identical (i.e. , segments 
do not overlap) . The innermost vrltie must 
be the first value of RS. Careful selection 
of the segment center locations is necessary 
to avoid introducing numerical inaccuracies 
associated with the use of finite vortex 
filaments. Smaller segments in regions of large 
bound circulation gradients (e.g*, at the 
blade tip) are required. Segment lengths 
less than 0.02R should normally be avoided. 

See Fig. 3 and the section* Special Input and 
Operating Instructions for a typical se^sat 
distribution. 
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LOCATIONS PROGRAM 
Symbol 

Deck Deck in 

F309 Yl4l» DecK F389 DESCRIPTION OF INPITT ITEM 

(Linked Version: prescribed internally as 
0.225, 0.375, 0.525, O.65, 0.75. O.85, 
0.925, 0.965, 0.99. First and/or second 
segments are Internally adjusted to conform 
to root cutout as defined by Location 133 
in Deck YlUl. ) 


39 STNEW Nujober of blade segments for the interpola- 

tion of the inflov solution and harmonic 
analysis of the inflov. If this it< . is 
zero the inflov at the solution segment 
centers (RS) are used for the harmonic 
analysis. li the sectional airfoil data 
option (Liication 207) is used these a**e the 
stations for vhich the sectional airfoil data 
are input. A maxin\u& of 25 is alloved. 
(Linked Version: Transferred from Deck YlHl 

as 15.) 

UO-6I1 RSNEtf Radial coordinates of blade segment centers 

of Interpolation calculations (nondimension- 
alized by R) . These segments are not lised in 
the circulation and inflov calculations. 

They are used only to interpolate the Inflov 
solutions at the original segments (RS) to 
these nev radial locations. Normally, RSNDf 
values are not input and the original RS 
se^aents are xised throughout the computation. 
A limit of 25 RSNEW values exists in the 
program. The innermost value must be the 
first value in RSNEW. (Linked Version: 
calculated in Deck YlUl from information In 
Locations 5-19 and transferred.) 
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DESCRIPTION OF INPUT ITEM 


VIMOM 


307 ALPHAS 


515 THET75 


136 THETAl 


Wake transport velocity, fps. Positive 
upflov. This velocity is only used if the 
undistorted wake model is to be internaJLly 
generated in the program. It is used to 
define the transport velocity, normal to 
the rotor disc. Normally, the momentum 
inflow velocity is used as defined by the 
following equations: 

I CtAB 

VIMOM = - i T 

where R, Op, , and are the tip speed, 
thrust coefficient, advance ratio and inflow 
ratio as conventionally defined in Ref. 7* 

(Linked Version: calculated in Deck YlUl and 
transferred. ) 

Rotor shaft angle, deg. Positive shaft 
angle corresponds to a nose->up tilt of the 
aircraft. Referenced to the forward fli^t 
velocity. 

Blade collective pitch at the 0.75R station, 
deg. This item is used in conjunction with 
THETAl and DTHETA to calculate the blade 
pitch distribution. 

Blade linear twist rate, deg. This input 
item is used to calculate the blade built-in 
pitch distribution for linearly twisted 
blades, and is normally negative to indicate 
decreased built-in pitch angles at the blade 
tip. The blade pitch, THETA, at each station, 

I, is calculated from the following equation: 

THETA(I) ■ THET75 ♦ THETAl (RS(l)-0.75)+DTHm(l) 


, ta 99 »rw 0 





Sikorsky PIrcraft 


REPORT NO. SER-50910 



LOCATIONS 

PROGRAM 



Deck 

Deck 

Symbol 

in 

DESCRIPTION OF INPUT 
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F38? 

YlUl* 

Deck F389 



69-83 


DTHETA 

Increment of blade pitch angle 

at each inflow 


solution blade station (RS), deg. For a 
nonlinear twisted blade, this item is used in 
place of or in conjunction with THETAl to define 
the built-in pitch angle at each station, I, 
due to blade twist (see above equation). The 
value for the innermost station must be the 
first value of DTHETA. (Linked Version: 
interpolated internally at radial stations for 
inflow solution (Locs. 24-38) from information 
in Locations l65 to 179 in Deck Yl4l.) 


84 

513 

AISP 

The first cosine harmonic coefficient of cyclic 
pitch, deg. Shaft axis system. Based on 
negative Fourier series. 

85 

514 

BISP 

The first sine heurmonic coefficient of cyclic 
pitch, deg. Shaft axis system. Based on 
negative Fourier series. 

86 


DELPSI 

Blade phase angle, deg. This item can be used 
to position the rotor in any desired starting 
location relative to the zero azimuth position. 
It is positive in the direction of rotation. 
Normally not input and blade 1 of the rotor is 
assumed to lie at ?ero deg for rotor position 1. 
(Linked Version: prescribed internally as 0.0} . 

87 


DELTA 

Bl'ide lag-angle, deg. Positive opposite to 
the direction of rotation. A constant lag 
angle is assumed. (Linked Version: 
transferred from Deck Yl4l as the mean lag 
angle computed therein.) 

88.96 



Not used. 

99 



Multiple case control card, see section entitled 


Multiple Case Control and Secondary Variable 
Input. (Linked Version: multiple case control 

cards must be provided.) 
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DESCRIPTION OF INPUT ITO 


100-120 


121-lkl 


1U2-162 


CURVE Airfoil lift curve slope vs. Mach numoer table, 
Cj^/rad. The first location is the nurnher of 
pairs of data. The data pairs follov 
(Mach No., C]^/rad) vith a minimum of 2. The 
data pairs are ordered in Mach number, lowest 
Mach number pair first. An interpolation 
routine is provided in the program which 
linearly interpolates on Mach number. Care 
must be taken to cover the Mach number range 
involved. An initial negative Mach number is 
required for the reverse flow region. These 
input values are used in the linearized 
circxilation solution. (Linked Version: not 

used in linked mode. This information is 
provided by lift curve slopes at each segment 
and azimuth as trcuisferred from Deck YlUl. 

See section on Sectional .^irfoil Characteristics 
Input.) 

ACURVE Angle of maximum lift coefficient vs. Heich 
number table. Uses the same input format as 
CURVE. These input vedues are used to define 
the stall limits in the linearized circulation 
solution. (Linked Version: not used in linked 

mode. This information is provided by angles 
of maximum lift coefficient at each segment and 
azimuth as transferred from Deck YlUl. See 
section on Sectional Airfoil Characteristics 
Input.) 

AOCURV Angle of zero lift vs. Mach number table, deg. 

For blades with cambered airfoil sections. Uses 
the same input format as CURVE. These input 
values are used in the linearized circulation 
solution. (Linked Version: not used in linked 

node. This information is provided by angles 
of zero lift at each section and azimuth as 
transferred from Deck Yl4l. See section on 
Sectional Airfoil CharflLCteristics Input.) 


3kS. 

PAOC 


RCV 0 


Sikorsky PIrcraft 


REPORT NO. SER-5O91O 


CMVMH umtvo 



eOHMAATlON 


LOCATIONS PROGRAM 
Symbol 

Deck Deck in 

F389 YlUl» Deck F389 DESCRIPTION OF INPUT ITEM 

and a 4.0 indicates a wake model from the tape 
output of another program. (Linked Version: 
prescribed internally as 1.0.) 

187 RUN Run number for wake geometry location on tai>e. 

A 1.0 is \ised for TYWAKE=1., 2., and 3. 

(Linked Version: prescribed internally as 1.0*) 

188 TSTEP Time step in tape file to start iislng vake 

geometry. Used only when TYWAKE = 4. 

(Linked Version; prescribed internally as 0.0.) 

189 REV Number of wake revolutions. Selection of this 

item is a compromise between accuracy and 
computer cost. The following values are typical: 
for hover (advance ratio, O) use at least 
0 revs; for 0— 0.05 use 6 revs; for 0.05 ~ 0.1 
use 4 revs; for 0.1 — 0.15 use 3 revs; for 
fL 2 0.15 use 2 revs. (Linked Version: 
prescribed internally as 5 0. Caution - it may 
be desirable to overwrite this value as 
described above.) 


190 ROLLUP The number of filaments to rollup in the tip 

vortex (at rollup angle TRUHC). If rollup is 
not desired, such as for a classical wake, use 
0.0 or 1.0. If rollup is desired, this value 
should be the number of tip vortex filaments 
outboard of the radial station at which the 
peak circulation occurs. This value is assumed 
constant (not varying with azimuth position), 
and thus a representative veilue must be 
selected. Since this value is not known a 
priori, it must be estimated, and if necessary 
iterated on. In hover the peak circulation 
normally occurs between 0.9 and O.90R* 
high speed forward flight it normally occurs 
between 0.7 and 0.85R. (Linked Version: 
prescribed internally as 5.0. Caution - It may 
be necessary to overwrite this value as 
described above.) 
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F389 Yllg* Deck F389 DESCRIPTION OF INPlfT ITEM 


191 AIJDTIP Option indicator for anhedral tip design. 

A O.t iriieater no droop, a 1.0 indicates a 
linear droop computed interneilly using DROOP, 
and a 2.0 Indicates input droop cooi-dinates 
(see Loc. 195). (Linked Version: prescribed 

internally as 0.0.) 

192 DROOP Linear droop aingle for anhedral tip, positive 

down, deg. (Linked Version: prescribed 

internally as 0.0.) 

193 RBEND Radial location for the start of the droop 

section, nondimensionalized by R. It must be 
positioned such that there are no more than 
5 inflow stations outboard of this location. 
(Linked Version: prescribed internally as 0.0.) 

19U RATIOV Ratio of the tip filament displacement velocity 

to VIMOM of the wake. If not input or a 0.0 la 
input the value is automatically set equal to 
1.0. (Linked Version: prescribed internally 

as 0.0). 

195-199 ZANDTP Axial coordinates of the drooped stations 

nondimensionalized by R. Referenced to unconed, 
undrooped blade. Dihedral: positive; 

Anhedral: negative. Maximum of 5 are allowed. 

This is used only if ANDTIP is set to 2.0. 

(Linked Version: prescribed internally as 0.0.) 

200 TRUNC Angle of rollup of the tip vortex filaments 

with respect to the blade, deg. Tip vortex 
filaments beyond this wake azimuth angle are 
all assigned coordinates equivalent to the 
outermost tip filament, thereby simulating 
the rollup of the tip vorticity into a concen- 
trated tip vortex. TRUNC/DPSI must be an integer* 
A typical value of TRUNC is 30. If no rollup is 
desired, set TRUNC « 0.0. (Linked Version: 
set equal to DPSI (Location U). Transferred 
from Deck YlUl). 
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F389 YlUl» Deck F389 DESCRIPT~ON OF INPUT ITEM 


201 


202 


203 


20U 


TRUNCI Truncation angle of vortex sheet filaments 
with respect to the blade, deg. Normally, 

TRUNCI * 0.0 which suppresses the truncation 
featxire. If a truncation angle is desired, 
TRUNCI = wake azimuth angle for truncation. 
(Linked Version: prescribed intenially as 

0 . 0 .) 

RCORE Tip vortex filament core radius nondimension^ 
all zed by the rotor radius. Generally a value 
of approximately 0.1 of the blade chord is used. 
Caution - RCORE must not exceed one-half of 
the length of the outer blade segment. (Linked 
Version: calculated internally as one-tenth the 

tip chord.) 

RCOREI Vortex filament core radius for filaments 
representing the inboard vertex sheet, 
nondimensionalized by the rotor radius. 

Caution - RCOREI must not exceed one-half 
the length of the smallest blade segment. 

(Linked Version: calculated internally as 

one-tenth the tip chord.) 

Not used. 


20$ TOL Tolerance control for the matrix solution of 

circulations. Normally, set TOL » 0.0 euad a 
tolerance of 0.000$ will be used (tolerance ■ 
allowable difference in suiamation of circulation 
solutions between last and previous iteration 
divided by summation for last iteration). 

Thus 0.0005 represents a 0.05% accuracy of 
convergence of the solution summation. If a 
tolerance value other than 0.0005 is desired, 
set TOL equal to the desired value. (Linked 
Version: prescribed internally as 0.0.) 
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F389 YlUl* Deck F389 


206 


207 


203 


209 


210 



211 
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WIOPT 


WOPT 


OPTPO 


PRINT 


PUNCH 


CPOPT 


Option indicator for W1 input. If a value 
of 0.0 is input W1 is computed from rigid 
blade equations using flapping harmonics and 
blade control parameters. A value of 1,0 
indicates that Wl's will be read in via cards. 
See section entitled, Noninduced Axial Velocity 
Component Input, Wl. (Linked Version: 
prescribed internally as 1.0. Input cards are 
replaced by data transferred from Deck YlUl), 

Option for sectional airfoil data input via 
cards. A 0.0. indicates no unsteady 
aerodynamics. A value of 1.0 indicates that 
sectional airfoil data input will be read from 
cards. See section entitled. Sectional Airfoil 
Characteristics Input. (linked Version: 
prescrioed internally as 1.0. Input cards are 
X <.-pla.ced by data transferred from Deck Yl4l. ) 

Geometric coefficients print option. A 0.0. 
indicates no print; a 1.0 calls for the print 
option. (Linked Version: prescribed Internally 

as 0.0.) 

The matrix solution print option. A 0.0 indi- 
cates no printout; a 1.0 calls for printout 
of the matrix coefficients and solution at 
each iteration, and a 2.0 prints only the 
solutions at each iteration. (Linked Version: 
prescribed internally as 0.0.) 

Card punch option for induced velocities. A 
0.0 indicates no punch cards, a 1.0 requests 
harmonics of induced velocity, a 2.0 requests 
Induced velocities at each point in the rotor 
disc, and a 3>0 requests both harmonics and 
point values. (Linked Version: prescribed 

internally as 0.0.) 

Make coordinate print option. A 0.0 indicates 
no print, a 1.0 calls for a printout of the wake 
coordinates. (Linked Version: prescribed 
.internally as 0.0.) 
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DESCRIPTION OF INPUT ITEM 

212 

CRdPT 

Print option for vake-element point 
combinations within the core radius. A 0.0 
input indicates no printout, while a 1.0 calls 
for the printout. (Linked Version: prescribed 

internally as 0.0.) 

213 

CPUNCH 

Option for punched cards of the rotor 
circulation solutions. A 0.0 indicates no 
punch, while a 1.0 calls for punched card 
output. (Linked Version: prescribed 

internally as 0.0.) 

214-219 


Not used. 

220 

OPLINK 

Option for running linked to other aeroelastie 


type programs. A 0.0 indicates normal 
unlinked operation. A 1.0 indicates normal 
linked operation, and the program vill 
automatically adjust for root cutout. If the 
root cutout is inhoard of the second inflow 
station, OPLINK should he set equal to 2.0. 

This vill bypass the root cutout adjustment vhen 
running linked, hut it also implies that the 
required inflow stations vill he input via the 
secondary loader input. (Linked Version: 
prescribed internally as 1.0.) 

221 DESUGP Option for intermediate printout, a 0.0 indicates 

no printout, a 1.0 requests printout. (Linked 
Version: prescribed internally as 0.0.) 
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Description of Input Item 


13-2U, 

25-36, 

37-U8, 

1*9-60. 


Starting at the innermost station, 5 
to a card for as many stations as are 
required. 

3 

• 


TAB 

e 


• 

• 

B 


• 

TAB 


H+1 

1-12 

13-2U, 

25-36. 

37-1*0. 

1*9-60 

TAB 

Angle of zero lift, starting at the 
Innermost station, 5 to a card for ais 
many stations as are required (deg.). 

B+2 

• 


TAB 

• 


• 

m 

n*H 


• 

• 

TAB 


N+Mfl 

1-12, 

13-21*, 

25-36. 

37-1*0, 

U9-6O 

TAB 

Slope of the lift curve, starting at 
innermost station, 5 to a card for as 
stations as are required (Cj^/rad.). 

B+M4-2 

• 


TAB 

• 


• 


• 

TAB 



Repeat the sequence starting vith card 1 for the next hlade position » until 
data for all rotor azimuth positions have been input. 
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Noninduced Axial Velocity Component Input, W1 


The ncnlnduced axial velocity component, Wl, m^ he either 
calculated internally in the program or read in via cards. The option 
control in the basic variable input is WIOFT. If calculated internally, 
blade flexibility is neglected and rigid, articulated blade equations 
are used to include the influence of rotor attitude, pitch angle and 
flapping motions, Wl values for each point In the rotor disc may be input 
on cards to include blade flexibility effects. This is particularly useful 
when the Wl values are calculated in another computer program. Appropriate 
values for eurticulated or nonart iculated rotors may be used. When 
including flexibility effects, the total noninduced axial velocity 
component, Wl, consists of axial velocity components due to both 
torsional and bending flexibility in addition to the rigid blade axlad. 
velocity. The Wl velocity is defined to be positive for upflov. 

The Wl values are input for the RSHEW stations and linearly interpolated 
in the program to the inflow stations, RS. The required input format is 
presented below. 

The data is input six items to a card using a Fortran 6E12.5 format. 
There is a set of cards for each radial interx>olated inflow solution 
station (RSHEW). The sets are ordered from the innermost to the outermost 
station. Each set of cards consists of the Wl values at all azimuthal 
positions for a given radial station. These values must start at the first 
actual blade azimuthal position (J » 1) , and are ordered azimuthally 
in the direction of rotation; therefore, care must be taken when blade 
phasing is used to insure the correct order of input. Stcurting with the 
first set, a card is filled with 6 items, then a new card is started; 
if needed. This is continued until all azimuthal positions for a given 
set are complete. For the next station this process is repeated. There 
are no title cards or blank cards used. For example, a case with 
1^ azimuthal positions (DPSI > 2U deg) and 9 blade stations (STNS *9) 
would have a total of 27 cards input for that rotor (3 cards per set x 
9 stations). 
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Tip Vortex Coordinate Overwrite Input 


In order to include a more realistic tip vortex geometry than the 
classiced \indistorted geometry, a tip vortex coordinate overwrite option 
is available. This option allows the program operator to input, via 
cairds, the tip filament geometry. This geometry can be obtained from 
experimental data or analytical results . With this option more realistic 
rear field effects can be obtained in cases where the wake is known to 
distort significantly from the classical model, and the geometry of the 
tip vortex is critical to obtaining accvirate inflow solutions. The required 
input format is presented below. 

The tip filament wake geometry is defined by inputting the vortex 
segment end point coordinates in the rotor wake reference system, (hub 
centered, parallel to the tip path plane and aligned with the free 
stream velocity) . Since the wake geometry is assumed to be periodic with 
blade spacing, the number of rotor positions is just 360/(B*D?SI) . In the 
program this quantity is called BJTOT. Thus, with the wake input for each 
blade at *0.1 of the rotor positions, the total wake geometry for all 
azimuthal positions is defined. The wake is input in sets of data where 
each set is associated with a rotor position (L). Each set is then 
composed of subsets of data corresponding to each blade at a given rotor 
position. The format is as follows: 

For each subset the first card describes the blade number (IB) and 
the nijmber of vortex segment endpoints to be input for that blade at that 
rotor position (NPM). These must be two integer, right adjusted numbers 
described by the Fortran format 215 (Columns 1-5, 6-10). On the following 
cards the x, y and z coordinates, in that order, for each endpoint are 
punched, • . 

The point input sequence is in the following order. The last point 
in the wake to be overwritten is input first, the next to last next, and 
so forth until the point corresponding to the blade segment boundary is 
input, this being the last point to be input for that filament. This allows 
the computer to retain sd.1 the classical wake points past the overwrite. 

Thus any portion of the tip filaments may be overwritten. The number of 
overwrite points for any filament must not exceed the total number 
associated with the corresponding classical wake. These items are input 
using a floating point format, eight spaces to an item, nine items to a 
card (3 segment end points). When this is completed the following blade 
subsets (ordered first to last) are punched in the same manner.* 

* The rotor blades su*e ordered as follows. Blade one lies at the first 

azimuthal position (including phasing) and the following blades (2 throu^ 
B) are defined to exist atimuthally behind the first blade. For examples 
with a 3-bladed rotor. If blade one is at zero azimuth, blade two is at 
2U0 deg and blade three is at 120 deg. 

ii£ 

2> .tv 0 







Sikorsky Pircraft 


REPORT NO. SER>50910 



To complete a set of data (one rotor position), this grouping of 
cards is proceeded by a single card describing the particular rotor 
position this set is for. A fixed point, right Justified format is used, 
Fortran format 110 (Columns 1-10 ). This format is repeated for each rotor 
position (ordered first to last) to complete the wake overwrite input. 
Figure 6 is a pictoral representation of this input sequence for one rotor. 
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Figure 6. Saiqple Cerd Input for Tip Vortex Geometry Overwrite 
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Sunvary List of Available Options 


In order to group together the available options for use in this 
program, the following list is presented. All of the items listed are 
controlled by input items in the basic vcuriable input. Listed are the 
available option controls and their locations. Details of their use wore 
described in the input Instructions. 


OPTIOH 


NAME OF 

LOCATIOH OPTION COlfTROL 


Interpolated induced velocities 

39 

STNSW 

Different types of wakes 

186 

TYVAKE 

Anhedral Tip 

191 

ANDTIP 

Wake tn’rcation 

200, 210 

TRUNC, TRUNCI 

Rollup of filaments into tip filamenc 

190 

ROLLUP 

Input of VI (nonlnduced inflow) 

206 

WIOPT 

Airfoil Characteristics 

207 

WOPT 

Geometric Coefficient print option 

208 

OPTPO 

Matrix Solution print option 

209 

PRINT 

Induced velocity punch option 

210 

PUNCH 

Wake coordinate print option 

211 

CPOPT 

Print option for blade points within 
vortex core 

212 

CBOPT 

Circulation punch option 

213 

CFURCH 

Program link option 

220 

OPLinC 

Intermediate print option 

221 

DEBUOP 
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DESCRIPTION OF PROGRAM OUTPUT 
Categories and Sequence of Program Output 


The categories of program output vhlch are printed In the foUovlng 
sequence are: 

1. Print of Basic Variable Input Data 

2. Intermediate Printout (optioned.) 

3. Wake Coordinate Printout (optional) 

U. Circulation Matrix Printout (optional) 

5 , Print of Circulation and Induced Velocity Solutions and Related 
Parameters 

6 . Print of Harmonic Coefficients of Induced Velocity 

In addition, punched card output of the induced velocity and circulation 
solutions is optional. A detailed description of the output items 
included in each of the above categories follows below. A sample computer 
printout of the output items is presented in a later section of this 
report. 

Print of Basic Vazlable Input Data 


This output is the same data which is input in the Basic Variable 
Input section via LOADER format. It is printed with a special output 
format to facilitate the program operator's ability to quickly change 
variable items without the necessity of using the program documentation 
once a basic understanding of the program is obtained. Descriptions of 
any of the items can be found in the input section of this report. The 
print format for each particzilar item is as follows. The first output 
quantity is the location of the item in the loader array. The second 
quantity is the variable name of the item as used in the program. A 
brief description of the Item and its units follows. The last quantity 
output Is the actual value of the item. 
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Intermediate Printout (Optional) 


When rxuining this program the user may or may not be interested in 
certain qu£jitities related to the solution of the particxaar case in hand. ] 

In order to keep printout to a minimum, the printout of this information 
has been made optional. The option control is DEBUCP in location 221 of 
the basic variable input. The first of these queuitities are the rotor 
blade coordinates in the cartesian coordinate system. These coordinates 
are listed in tabular form, horizontally as a function of radial station (RS) 
and vertically as a function of azimuth position (PSI). 

The following items are aerodvramic in nature and are local blade 
element quantities. They use the same output format as noted above, and 
are output in the following order. 

1. In-plane velocity component table, UT (fps). 

2. Mach nxmber table, MACH. 

3. Airfoil lift curve slope-table AT (Ci/rad). 

U. Angle of attack for maximum lift table, ALMAX (deg). 

5 . Angle of attack for zero lift table, AIX) (deg). 

6. Total axial noninduced velocity component table, W1 (fps). 

7 . The product of the lift curve slope and the inplane velocity 
component table, ACITTAB (ft2/sec). 

8. The maximum circulation table, based on stall considerations, 

CIRCM (ft2/sec). 
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VeJce Coordinate Printout (Optional) 


If desired, the coordinates of the vortex filament segment endpoints 
are printed for each filament string of each blade at each rotor position. 
The option control is CFOFT in location 211 of the basic variable input. 
The coordinates listed are nondimens ionali zed by the rotor radius. For 
any one filament, the first coordinate set is the last segment endpoint 
of the string (i.e., the oldest segment). The last value listed in the 
filament segment endpoint on the blade. All points in between are 
naturedly sequenced in order (i.e., oldest to newest). Tbe coordinate 
print is frequently omitted because of the printout page requirements. 


Circulation Matrix Printout (Optioned) 


The formulation of the circulation matrix coefficients was 
previously presented in the description associated with Equations 1 
throxigh 17. In summary, the geometric coefficients are defined as the 
geometric functions in the Biot-Savai*t relation used to calculate the 
velocities induced by the wake and bound vortex segments at the blade. 

Since the circxilation strengths of all wake segments within a given trailing 
vortex filament are assumed periodic with azimuth position for the forward 
fli^t condition, the geometric coefficients of a group of segments within 
a filament for the sane azimuthal position at the time of generation may 
be summed over all revolutions to establish a representative geometric 
coefficient for the group (GCTBOW). To relate these coefficients to the 
respective bound circulations, a difference is computed between the wake 
coefficients of the trailing filaments to obtain the geometric coefficients 
corresponding to the azimuthally veurying bound circulations as related to 
the wake structure. To these values are added the effects of the bound 
circulations (GCB). The res'ilting quantities, GCM/^T are used to form the 
coefficients of the circulation matrix as shown on p. in Equation 16. 
Within the matrix solving subroutines (MILSOP or MILS0P2), the matrix is 
normalized by dividing each of these values by the diagonal element. The 
option control for printing these circulation matrix coefficients is PRIBT 
(location 209) in the basic variable input. An option is edso provided by 
PRINT to output the intermediate circulation solutions at each iteration. 
These two options are generally not used because of the lengthy print 
requirements. If the circulation matrix coefficients are output, the 
matrix row and column is printed for each line of data for reference. 

The intermediate cirexilatlon solutions are output with it's position 
in the column vector. If desiied, the GCMAT quantities may be printed. 

The option control in the basic variable input is OPTPO (location 208). 
When this option is used thw output is in the following form. At each 
inflow point on the rotor disk the geometric coefficients are listed 
(10 to a line). They are ordered according to radial station for each 
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azimuth with no indication where the break from one azimuth position to 
another is made. This printout is not in as convenient an output form 
as the matrix coefficients and is not generally used. 

Print of Circxilation and Induced Velocity 
Solutions and Related Parameters 


The blade bound circulation and induced velocity solutions are 
printed for each radial station and azimuth position. These output items 
are not optional since they are the desired results of the program operation. 
Following the circulation and induced velocity printout, the blade angle of 
attack distribution based on Equation is printed. A stall indicator 
table is printed next. In this table, a 0.0 or a 99-99 is printed for each 
radial station and azimuth position. A 0.0 indicates no st«J.l. A 99«99 
indicates that the stall limit (ACL'RVE, see p. ) was exceeded and that the 
circulation at that station was limited to a maximum value in the 
circulation solution. If used, the interjxslated induced velocity 
distribution is then printed. 

Print of Harmonic Coefficients of 
Induced Velocity 


The induced velocity solutions are harmonically analyzed and the 
harmonic coefficients, based on a positive Fourier series, are printed. 

The harmonics are printed for the inflow stations (RS) unless interpolated 
radial stations (RSNEW) have been requested. 

T\inched Cards of Induced Velocity 
and Circulation Solutions 


An option for punched cards of the Induced velocity solutions at 
each blade radial station and azimuth position is available and controlled 
by PUNCH (location 210 ) in the basic vsuriable input. An option for punched 
cards of the harmonic coefficients of induced velocity at each blade radiad 
station is also controlled by PUIICH. An option for punched cards of the 
circulation solutions is controlled by CPUNCH (location 213). 



X 


I 


a 


SA 29 REV 0 



51c 

PAGE 


4 



I 


Sikorsky ?^lrcraft 



WVIMN UNmo oomaation 

fi. 


REPORT NO. SER-50910 


SPECIAL INPUT AND OPERATING INSTRUCI'IONS 


This program will allow the use of a wide range of rotor blade 
geometries. Variable twist blades (both linear and nonlinear) Cfui be used. 
Variable chord distributions can be used. No restriction on the ntonber 
of blades is made. The program assxnnes solution periodicELLly with blade 
spacing. The choice of DPSI and the number of blades is interdependent 
since DPSI must be such that the number of rotor positions between blades 
and the number of blades are compatible (360/B/DPSI = integer). The 
flapping hinge offset must be inboard of the innermost vortex sheet 
filament . 

Concerning th'^ choice of radial stations for the inflow solutions* 
careful selection cf stations is necessary. Since the inflow stations 
also define the ti ailing vortex filament spacing, attention must be given 
during their selection to the wake model that resxilts. Generally the tip 
region should contain a larger number of filaments per iJiit length than the 
inboard region. This allows a better modeling of the concentration of 
trailing vorticity in the region of the tip vortex and also gives better 
resolution of the bound circulation nesur the tip (this is needed because of 
the severe loading gradients observed in this region) . Standard values which 
are generally used at UARL for a 9 station case can be obtained from the 
input description of RS in the basic veiriable input section of this appendix. 
The inboard values must be adjusted for different root cutouts of various 
blades. If an anhedral tip is being used* a maximum of five Inflow solution 
stations is allowed outboard of the start of the anhedral tip section. 


When modeling wake rollup in forweurd flight* ceureful selection of 
the number of filaments that rollup into a concentrated tip vortex is also 
necessary. The selection is based on taking the average of the radial 
position of peak circulation on the disk and choosing the number of 
filaments for rollup to be the total of those which are outboard of this 
average radial position. Since the locations of the peak circulations are 
not known beforehand an Iteration procedure may be necessary. Generally 
for the stations mentioned above, the 5 outboard filaments are "on the 
average" outboard of the average radial location of peak circulation. 


SAMIIIE!V« 


The blade motion and control parameters should be chosen to be 
compatible with the wake model used, since these parameters determine to 
a large degree the rotor thrust, blade loading and other parameters which 
are directly related to the wake geometry. Generally a separate program 
is used to compute these parameters for a desired operating condition. 
With these parameters, other appropriate input, and a wake model, the 
resulting inflow distribution la computed with this program. This inflow 
is then used in the program which generated the control parameters values 
and rerun. If the new control parameters are reasonably close to the 
required values the wake model and control parameters cure compatible. If 
not, more iterations may be necessary. 
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Operation of this program is gene, ally straightforward once a basic 
understanding of the program input is c btained. However, there is one area 
where problems may occur. For cases in which vortex elements come 
within extreme close proximity to blade.>, large xuireollstic geometric 
coefficients can result. These large values can lead to large off-diagonal 
terms in the circulation matrix which results in convergence problems in 
the matrix solution or locally unrealistic solutions. This may generally 
be avoided by proper selection of the core radius assigned to the vortex 
filaments. (The velocity induced by vortex element is set to zero for 
blade points wihtin the vortex core radius). Experimental data indicate 
tip vortex core radii of approximately one-tenth of the blade chord for 
helicopter rotors. However, the proper core radius for the inboard 
filaments representing the vortex sheets is not well defined at this time. 
\?alues of one-tenth to three-tenths of the blade chord have generally 
been used. However, some unrealistic solutions associated with close 
proximity have resxilted in certain instances with the former value. 
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